KAGUYA LISM MAP DTM_MAP import into ISIS3
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?
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
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
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
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
#2 Updated by Tammy Becker over 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:
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).
This can be clipped and reprojected in MAP2 (with an account) if you can't download the 21GB ISIS3 cube!
Hope this helps!