Jump to content
PDS Geosciences Node Community

June Wang

Geo Staff
 View Profile  See their activity

  • Posts

    616

  • Joined

  • Last visited

 Content Type 

Posts posted by June Wang

    Convert 64-bit MRO RSS Digital Map in PDS IMG format to 32-bit GeoTiff with GDAL

    in For data users

    Posted

    We were contacted by a user who are looking at the tectonics and gravity of the Tharsis Rise region of Mars. She mentioned that she has “run into a problem analyzing some PDS data that I can't find a way around. After many days of trial and error (which was not aided by the fact that I have very little computer science experience), I was hoping you might be able to help. 

    I'm trying to view and analyze the GMM-3 Mars gravity models, which are .IMG files (PDS listing). I don't have access to ENVI, so I've been trying alternatives, and have had the most success with the USGS's ISIS software. Now I'm stuck getting the IMG file into an image cube format that ISIS can display and analyze using qview. The pds2isis and raw2isis applications won't work because they can only process up to 32-bit images, and these are 64-bit (label file).  

    So far I have tried, and not had any luck with,

    •  Converting the 64-bit IMG files to 32-bit files

    • Finding lower-quality versions of the data

    • Finding alternatives to ISIS that could process or display the IMG files

    My main challenge is that I need to view the data with its original unit values in milligals. Programs like ImageJ and Photoshop can open the raw PDS image files, but it will lose the units of the raw data, replacing them with a standard grayscale with pixel values of 0-255 that are disassociated from the original data values in milligals. I had the same problem with the software NASAView. 

    If it's relevant, what I want to do with the data is create maps and profile plots across a transect for gravity and elevation data that have the actual unit values so I will be able to include a map scale and real quantitative values on the plots.”

    We suggested below to use GDAL for working with MRO RSS Digital Map conversion.

     “GDAL would help you to translate the img data to another format, but the original PDS label defined “MAP_SCALE  = "N/A" “ in the *.lbl file. GDAL probably didn’t like this and didn’t work for me in my first try. I manually changed that string to ‘’ MAP_SCALE                    = 3704.461337 <m/pixel>” as shown in the attached label, and run below gdal_translate commands. GDAL could read the PDS IMG data and output a Geotiff. The first command output a 64bit Geotiff, and the second command output a 32 bit Geotiff. You can read the page (https://gdal.org/programs/gdal_translate.html) to adjust the parameters for your application.

    gdal_translate -of GTiff D:\test\data\mro\ggmro_120_anom_100.lbl D:\test\data\mro\ggmro_120_anom_100.tif

    gdal_translate -ot Float32 -of GTiff D:\test\data\mro\ggmro_120_anom_100.lbl D:\test\data\mro\ggmro_120_anom_100_32.tif

    GDAL can be downloaded from https://gdal.org/download.html . You can also find some tutorial on the web. I usually download GDAL for my windows computer from GISInternals (https://www.gisinternals.com/release.php). If you have a windows operating system, you can try the new version of “gdal-302-1928-x64-core.msi”near the middle of the page, https://www.gisinternals.com/query.html?content=filelist&file=release-1928-x64-gdal-3-2-1-mapserver-7-6-2.zip . Double click the *.msi, the installation window will pop-up and it is straightforward.”

    Georeferencing PDS images in ENVI

    in For data users

    Posted

    Yes, ArcGIS has better projection engine defined. And ENVI used ArcGIS projection engine in its application. The steps 3-9 you mentioned are actually based on ENVI classic interface. That blog was written years ago. After that there are many new release of ENVI, and its current version is 5.6. I don’t know which version you are using now, probably something 5.4+. If you are using a newer version of ENVI interface, one way to add the spatial reference is to find “Edit ENVI Header” under “Raster Management” in the toolbox as shown below in the right. You can add spatial reference from there. Instead of adding something for the moon, you will need to add the projection for Enceladus for this data product.

    image.png.e63e2b3418e19e9522b5593b13cd79d5.png

     

    When I read the label file, I noticed this data product is a CASSINI EDR data. It is archived at PDS Imaging Node under volume https://pds-imaging.jpl.nasa.gov/data/cassini/cassini_orbiter/coiss_2042/ . I would not suggest to use this EDR directly in your crater counting.

    If you read the “aareadme.txt” file in this volume, it said “The ISS team is not systematically performing calibration processing on the ISS EDRs.” However, their 'CALIBRATION' volumes (COISS_0xxx) “include all the necessary tools and documentation to support future users in performing their own calibration processing of the ISS EDR image datasets.” And the section “3.  Data Calibration“ of “aareadme.txt” has mentioned “The Cassini ISS archive volumes contain an extras directory that is included only on the eleventh Calibration volume. This directory contains the source code for the Cassini Imaging Science Subsystem Calibration (CISSCAL) software. This software, developed by the Cassini Imaging team, allows the user to radiometrically and geometrically process the EDR-level images into higher level calibrated images.”

    You can also read the SIS file (https://pds-imaging.jpl.nasa.gov/data/cassini/cassini_orbiter/coiss_2042/document/edrsis.pdf) for more information. If you contact with the people at PDS Imaging Node, or people from the CASSINI team, they may provide expert help for you on this topic than I.

    Thanks,

    June

    ODE - Release of PDS4 Migrated Lunar Prospector Data Archives at the Geosciences Node Loaded into ODE

    in Announcements

    Posted

    February 3, 2021
    Release of PDS4 migrated Lunar Prospector data archives at the Geosciences Node are loaded into ODE. This includes calibrated GRS and NS data, GRS elemental abundance tables, and RSS spherical harmonic gravity models archived in the below four bundles.
        lp_grsns_calibrated
        lp_grs_derived
        lp_rs_goddard_gravity_derived
        lp_rs_jpl_gravity_derived

    Please see below links for more information.
    https://pds-geosciences.wustl.edu/missions/lunarp/index.htm
    ODE Moon Holdings: https://ode.rsl.wustl.edu/odeholdings/Moon_holdings.html

     

    ODE - New MEX MARSIS EDR and Subsurface RDR Data Loaded into ODE

    in Announcements

    Posted

    January 21st, 2021 - New Extended Mission 7 MEX MARSIS EDR and Subsurface RDR Data loaded into ODE.
    - Updated MEX MARSIS EDR data products released through December 31st, 2019 (Orbit 1- 20228)
    - Updated MEX MARSIS Subsurface RDR data products released through December 29th, 2019 (Orbit 1885-20221)
    See https://ode.rsl.wustl.edu/odeholdings/Mars_holdings.html for a full list of data holdings in ODE.

     

    ODE - Mars Odyssey Release 74 and Re-release of 50 &51 Data Loaded into ODE

    in Announcements

    Posted

    January 7th, 2021 - Mars Odyssey Release 74 and Re-release of 50 &51 Data Updates
    Loaded NS and HEND raw and derived data into ODE through June 30th, 2020
    Updated THEMIS standard data products (IREDR, IRRDR, IRBTR, VISEDR, VISRDR and VISABR data) from February 19th, 2002 to May 31st, 2020
    Updated THEMIS geometry data products (IRGEO2, IRPBT2, VGEO2, and VISALB data) from September 28th, 2013 to May 31st, 2020
    Updated re-release volumes for THEMIS VIS calibrated and derived products from release 50 and 51

    See more information at https://ode.rsl.wustl.edu/odeholdings/Mars_holdings.html and http://static.mars.asu.edu/pds/errata.txt

     

    ODE - LRO Release 44 Loaded into ODE

    in Announcements

    Posted

    December 15th, 2020 - LRO Release 44 has been loaded into ODE.
    Updated LAMP EDR and RDR products from July 6th, 2009 to September 14th, 2020
    Updated LEND EDR_SCI, EDR_HK, RDR_RSCI, and RDR_CHK products from June 18th, 2009 to September 14th, 2020
    Updated LEND RDR_DLD and RDR_DLX products from June 27th, 2009 to June 14th, 2020
    Updated Diviner EDR and RDR products from July 5th, 2009 to September 16th, 2020
    Updated LROC EDR, CDR and RDR products from June 30th, 2009 to September 15th, 2020
    Loaded Mini-RF new Version 3 Bistatic radar data products (EDR, RDR, and DDR) from January 28th, 2019 to March 31st, 2020
    Updated LOLA EDR products from June 18th, 2009 to October 15th, 2020
    Updated LOLA RDR products from July 13th, 2009 to October 15th, 2020
        
    Updated LOLA RDR query tool database from July 13th, 2009 to October 15th, 2020
    In Progress Diviner RDR query tool database from July 5th, 2009 to September 16th, 2020

    Please see ODE Lunar Holdings - https://ode.rsl.wustl.edu/odeholdings/Moon_holdings.html

     

    ODE - MRO Release 55 Loaded into ODE

    in Announcements

    Posted

    December 3rd, 2020 - MRO Release 55 Loaded into ODE.
    -    Updated MRO CRISM EDR and CDR products from September 27th, 2006 to August 9th, 2020
    -    Updated MRO CRISM DDR products from September 27th, 2006 to August 8th, 2020
    -    Updated MRO CRISM TRDR products from September 27th, 2006 to August 8th, 2020
    -    Updated MRO raw gravity data products released through November 10th, 2020
    -    Updated MRO SHARAD EDR products from ASI through May 9th, 2020 (orbits through 64623)
    -    Updated SHARAD Radargram data products released through May 7th, 2020 (orbits through 64597)
    -    Updated MRO MCS EDR, RDR and DDR products from September 15th, 2006 to July 31st, 2020 (orbits through 65689)
    -    Updated MRO CTX EDR products from August 30th, 2006 to May 31st, 2020
    -    Updated MRO HiRISE EDR, RDR, DTM and anaglyph products released through November 3rd, 2020 (Orbit through 66899)
    Please see ODE Mars Holdings – https://ode.rsl.wustl.edu/odeholdings/Mars_holdings.html

    Note: SHARAD RDRs from ASI have not been delivered yet.

     

    ODE - Release 15 of MRO CRISM MTRDR, TER data is loaded into ODE, along with new TRDR browse and extras data products

    in Announcements

    Posted

    November 25th, 2020
    Release 15 of MRO CRISM MTRDR (Map-projected Targeted Reduced Data Record), TER (Targeted Empirical Record), and TRDR (Targeted Reduced Data Record) browse and extras data has been loaded into ODE. This release fills in data from August 9 through November 8, 2008.
    See https://ode.rsl.wustl.edu/odeholdings/Mars_holdings.html for more information.

     

    ODE - MRO Release 54 Loaded into ODE

    in Announcements

    Posted

    October 19th, 2020 - SHARAD EDR and RDR data from ASI for MRO Release 54 loaded into ODE, also loaded missing SHARAD RDRs from previous MRO Release 24.
    - Loaded delayed SHARAD EDRs and RDRs from release 54 (Data coverage: November 10, 2019, through February 8, 2020)
    - Loaded SHARAD RDR (ASI) missing data from release 24 (Data coverage: June 3 through August 3, 2012)

     

    Discrepancy in Line Projection offset value in PDS label for LRO NAC mosaic file

    in For data providers

    Posted

    One more finding is if using the function of “y = (1 - L0 - Line) * Scale” in line 128 of ‘DSMAP.CAT’ file, y value will always be negative with current label value input, which doesn’t match with the real y value in the data. It looks to me this equation could be changed to “y = ( L0 – Line-1) * Scale”.

    You can try LINE_PROJECTION_OFFSET = 4476. This value is just my guess, may not be correct.

    I have contacted with the data provider. I will post the final solution here when I hear anything back from the data provider’s.


    Details needed for LROC NAC mosaic .IMG file

    in For data users

    Posted

    Hi Ayushman,

    Two things here. In your code, you looped around all the pixels. Please don’t include the nodata area in your calculation.

    The LINE_PROJECTION_OFFSET value in the label doesn’t match with the value of CENTER_LATITUDE. Need to check with the data provider to see what the value used in their processing pipelines.

    Thanks,

    June

     

    Details needed for LROC NAC mosaic .IMG file

    in For data users

    Posted

    Hi Ayushman,  

    Attached is the screen capture of the LROC RDR (orange color, NAC_ROI_VIKRAM__LOC_P708S0237_5M.TIF) overlain on a WAC mosaic near the South Pole area together with its spatial reference information. The feature matches well. If everything is correct in the book, it looks to me equation 20-15 is a better fit for this case. Have you checked other parameters or steps in your calculation? If you can send me your code and your original value input for the calculation, I can take a further look. Thanks,

    June

    Capture_20201009.JPG

    Capture_20201009_spatial.JPG

    Details needed for LROC NAC mosaic .IMG file

    in For data users

    Posted

    Hi Ayushman,

    Question 1-> It looks to me C = 2* ARCTAN[P/(2*Rp)] is correct, when compared this with the equation (21-15) in book page 159 of the book “Snyder, J. P. (1987). Map Projections: A Working Manual. U.S. Geological Survey Professional Paper 1395. Washington, DC: United States Government Printing Office.” Please let me know if this is not correct to your test.

    Question 2-> Please use the value in the attached label. The sentence written in the DSMAP.CAT is more like a general case.

    Thanks,

    June

    Details needed for LROC NAC mosaic .IMG file

    in For data users

    Posted

    Hi Ayushman,

    To get the lat&lon for each pixel of the LROC map projected RDRs, please read the ‘DSMAP.CAT’ at https://pds.lroc.asu.edu/data/LRO-L-LROC-5-RDR-V1.0/LROLRC_2001/CATALOG/, which gives the equations and steps on how to convert image array coordinates (Sample, Line) to map coordinates (x, y), and then from map coordinates (x, y) to planetocentric coordinates (Lat, Lon) for different kind of projections. All the parameters needed for the calculation could be find in the attached label of this product. Please let me know if you need more information. Thanks,

    June

     

     

     

    Location coordinates needed for LRO NAC CDR images

    in For data users

    Posted

    Hi Ayushman,
    Sorry, these are LROC NAC CDR data products. They are not map projected, but you can still use Arcmap and Gdal to read the data. I confused them with the LROC RDR data in your other post at https://geoweb.rsl.wustl.edu/community/index.php?/topic/2804-details-needed-for-lroc-nac-mosaic-img-file/&tab=comments#comment-4594. I will keep a post there on how to get the lat&lon for each pixel of the LROC map projected RDRs.

    In the index table under INDEX directory of the LROC CDR volumes, e.g., https://pds.lroc.asu.edu/data/LRO-L-LROC-3-CDR-V1.0/LROLRC_10**/INDEX/, you will find the lat&lon coordinates of the center point and 4 corner points. If you want to get the lat&lon for each pixel in the image, you will probably need to try the NAIF SPICE tools https://naif.jpl.nasa.gov/naif/toolkit.html or USGS ISIS https://isis.astrogeology.usgs.gov/index.html to see if they have function for this.  

    Thanks,

    June

×
×
  • Create New...