Project

General

Profile

Question #4377

KAGUYA LISM MAP DTM_MAP import into ISIS3

Added by francesco zucca over 1 year ago. Updated 6 months ago.

Status:
Resolved
Priority:
Normal
Assignee:
Category:
Applications
Target version:
-
Test Reviewer:

Description

Hi at all,

I'm in trouble with the importing of KAGUYA LISM DTM_MAP (3x3 degrees) [both M or S] level2 map projected data, available from SELENE Data Archive and I probably I need some hints.
If I use pds2isis to import the dtm file [ pds2isis image=/Kaguja/DTM_M/data/test/DTM_MAP_02_N42E033N39E036SC.dtm to=/Kaguja/DTM_M/data/test/DTM_MAP_02_N42E033N39E036SC.cub] I get an **ERROR* No target name in labels***** and the output cub file is missing all the parameters (projection, file and data type..) I can get from catlab as at the end of message. As you can see there is not a TARGET NAME label, indeed.

Any idea about? Other command or sequence can help?

TIA
regards
fz

Output from catlab

PDS_VERSION_ID = PDS3
RECORD_TYPE = UNDEFINED
FILE_NAME = DTM_MAP_02_N42E033N39E036SC.dtm
PRODUCT_ID = DTM_MAP_02_N42E033N39E036SC
DATA_FORMAT = PDS
IMAGE = 3161
SOFTWARE_NAME = "LISM DTM generation software"
SOFTWARE_VERSION = 1.2
PROCESS_VERSION_ID = MAP
PRODUCT_CREATION_TIME = 2009-11-17T10:10:54Z
PRODUCER_ID = LISM
PRODUCT_SET_ID = DTM_MAP
PRODUCT_VERSION_ID = 02
MISSION_NAME = SELENE
SPACECRAFT_NAME = SELENE-M
DATA_SET_ID = DTM_MAP
INSTRUMENT_NAME = Terrain_Camera
INSTRUMENT_ID = TC
UPPER_LEFT_LATITUDE = 42.000000
UPPER_LEFT_LONGITUDE = 33.000000
UPPER_RIGHT_LATITUDE = 42.000000
UPPER_RIGHT_LONGITUDE = 35.999756
LOWER_LEFT_LATITUDE = 39.000244
LOWER_LEFT_LONGITUDE = 33.000000
LOWER_RIGHT_LATITUDE = 39.000244
LOWER_RIGHT_LONGITUDE = 35.999756
IMAGE_CENTER_LATITUDE = 40.500122
IMAGE_CENTER_LONGITUDE = 34.499878

Object = IMAGE_MAP_PROJECTION
MAP_PROJECTION_TYPE = "Simple Cylindrical"
COORDINATE_SYSTEM_TYPE = "BODY-FIXED ROTATING"
COORDINATE_SYSTEM_NAME = PLANETOCENTRIC
A_AXIS_RADIUS = 1737.400
B_AXIS_RADIUS = 1737.400
C_AXIS_RADIUS = 1737.400
FIRST_STANDARD_PARALLEL = N/A
SECOND_STANDARD_PARALLEL = N/A
POSITIVE_LONGITUDE_DIRECTION = EAST
CENTER_LATITUDE = 0.000000
CENTER_LONGITUDE = 0.000000
REFERENCE_LATITUDE = N/A
REFERENCE_LONGITUDE = N/A
LINE_FIRST_PIXEL = 1
LINE_LAST_PIXEL = 12288
SAMPLE_FIRST_PIXEL = 1
SAMPLE_LAST_PIXEL = 12288
MAP_PROJECTION_ROTATION = 0.0
MAP_RESOLUTION = 4096.000000
MAP_SCALE = 0.007403
MAXIMUM_LATITUDE = 42.000000
MINIMUM_LATITUDE = 39.000244
EASTERNMOST_LONGITUDE = 35.999756
WESTERNMOST_LONGITUDE = 33.000000
LINE_PROJECTION_OFFSET = 172032.000000
SAMPLE_PROJECTION_OFFSET = 135168.000000
RESAMPLING_METHOD = Bi-linear
End_Object

Object = PROCESSING_PARAMETERS
PARAMETER_SET_NAME = for_SCtoSCMAP
HORIZONTAL_TRANSFORM_METHOD = PARALLEL
VERTICAL_TRANSFORM_METHOD = TREND
MOSAIC_PRIORITY = (DATE_NEW, N/A)
SMOOTHING_WIDTH = 50
End_Object

Object = IMAGE
BANDS = 1
BAND_STORAGE_TYPE = BAND_SEQUENTIAL
BAND_NAME = N/A
LINES = 12288
LINE_SAMPLES = 12288
SAMPLE_TYPE = MSB_INTEGER
SAMPLE_BITS = 16
IMAGE_VALUE_TYPE = ELEVATION
SAMPLE_BIT_MASK = 2#1111111111111111#
OFFSET = 0.000000
SCALING_FACTOR = 1.000000
STRETCHED_FLAG = FALSE
VALID_MINIMUM = -9989
VALID_MAXIMUM = 32766
DUMMY = -9999
MINIMUM = -3925
MAXIMUM = -846
AVERAGE = -1895.941439
STDEV = 330.754678
MODE_PIXEL = -2030
End_Object

Object = QUALITY_INFO
QA_PERCENT_GOOD_PIXEL = 99.981986
QA_PERCENT_DUMMY_PIXEL = 0.018014
QA_PERCENT_BAD_PIXEL = 0.000000
QA_PERCENT_INTERPOLATED_PIXEL = 0.000000
QA_PERCENT_SHADOW_PIXEL = 0.000000
BAD_PIXEL_THRESHOLD_CORRELATION = 0.200000
BAD_PIXEL_THRESHOLD_SLOPE = 45.000000
End_Object
End

History

#1 Updated by Tammy Becker over 1 year ago

  • Status changed from New to Acknowledged

#2 Updated by Tammy Becker about 1 year ago

  • Assignee set to Tammy Becker

In order to ingest the dtm into ISIS3, the input label must contain the Target Name.

Trent Hare has offered the following suggestion:
Create a backup copy of the original dtm file. Replace an "unneeded" input keyword label with the TARGET_NAME.

Run the following PERL one-liner to swap RECORD_TYPE with TARGET_NAME, make sure to swap out the exact same number of characters (including blank spaces).

perl -pi -e 's/RECORD_TYPE = "UNDEFINED"/TARGET_NAME = "MOON" /g' DTM_MAP_02_N00E000S03E003SC.dtm

Now you should be able to run:

pds2isis from=DTM_MAP_02_N00E000S03E003SC.dtm to=test.cub

potential issues with created ISIS3 label.
(1) The longitude domain might be 180 instead of the default 360 as applied.
(2) pds2isis might introduced a 1 pixel offset in sample and line--this is a known issue. This potential ~7m error in X,Y (in the map projected Cartesian space) will show up in the UpperLeftCornerX and UpperLeftCornerY parameters in the ISIS3 label. For background information on this issue refer to:

https://isis.astrogeology.usgs.gov/fixit/boards/6/topics/8

In addition, Trent suggests to check out the version of the JAXA DTMs that have been corrected to LOLA. A collaboration by JAXA and Goddard. It helps with vertical and horizontal shifts in the original data. And while 59m/p sounds very sub-sampled from ~7m, it is really more realistic. The original DTM might be good at around 20 to ~30m but ~59m/p was a good compromise by Goddard (and helped with the registration and beat down errors).

http://astrogeology.usgs.gov/search/map/Moon/LRO/LOLA/Lunar_LRO_LrocKaguya_DEMmerge_60N60S_512ppd

This can be clipped and reprojected in MAP2 (with an account) if you can't download the 21GB ISIS3 cube!

Hope this helps!

#4 Updated by Tammy Becker 6 months ago

  • Status changed from Acknowledged to Resolved

Also available in: Atom PDF