ALOS PRISM CEOS Terrain Correction problems
Submitted by dyp on Tue, 10/30/2007 - 09:46
I am trying to do a terrain correction for ALOS PRISM granule in CEOS format. It is 1B2R Level. I have DEM from SRTM resampled to 2.5m. ASF version 3.2.9 as 3.1.16 says that it does not support CEOS terrain correction. When I save DEM to float GeoTIFF I get a crash: Exception: STATUS_ACCESS_VIOLATION at eip=610D8F8A eax=10027C58 ebx=0054CB56 ecx=00000000 edx=0054CB56 esi=00229544 edi=10029878 ebp=00229418 esp=00229414 program=C:\Program Files\ASF_Tools\Convert\asf_convert_gui.exe, pid 3132, thread main cs=001B ds=0023 es=0023 fs=003B gs=0000 ss=0023 Stack trace: Frame Function Args 00229418 610D8F8A (00000000, 0054CB56, 10028D88, 3FEFFCB9) 00229F68 00430740 (013649A0, 0022A3C0, 0022A410, 0022A518) 0022A538 004325C2 (013649A0, 0022AB60, 005436DE, 312D5755) 0022A7E8 00426F8F (00000001, 00000000, 00000000, 00000000) 0022AC78 0041A455 (1002A0B0, 00000001, 00000000, 00000000) 0022B838 0041BB95 (00000000, 1002A7C0, 00000001, 1002A3D0) 0022DA88 00402737 (01ADD734, 00000000, 00000003, 0022DAD4) 0022DB08 00405AA0 (017E7B50, 00000000, 0022DB48, 008C3806) 0022DB38 00AD85C2 (019CA4D0, 00000000, 00000001, 0022DD10) 0022DB78 00AC3B62 (019CA4D0, 00000000, 00000001, 0022DD10) 0022DC28 00AD7EF7 (017CDA98, 00000000, 017E7B50, 00000000) 0022DEA8 00AD7496 (017E7B50, 0000006D, 00000000, 0022DEDC) 0022DEC8 00AD7708 (017E7B50, 0000006D, 00000000, 00AE60FE) 0022DEF8 00BA77F7 (017E7B50, 017CD7F0, 0022DF38, 00AD85C2) 0022DF08 00BA8AB6 (017E7B50, 017916D8, 0022DF38, 008CF1FB) 0022DF38 00AD85C2 (017CD7F0, 00000000, 00000001, 0022E140) End of stack trace (more stack frames may be present) I have saved DEM as ENVI IMG file and now I have completely different error message: Starting terrain correction pre-processing. Reading SAR metadata from: C:\maps\inbox\space\GeoNadir\ALPSMW086962520-O1B2R_UW\
ALPSMW086962520-O1B2R_UW-1B2G-31-Oct-2007_00-14-26/import
** Error: ********
Invalid metadata for: C:\maps\inbox\space\GeoNadir\ALPSMW086962520-O1B2R_UW\ALPSMW086962520-
O1B2R_UW-1B2G-31-Oct-2007_00-14-26/import No SAR block found.
** End of error **
I can not understand what is the problem. Is it possible to make a terrain correction for ALOS PRISM data using ASF? Finally I want to have 1B2G Level granule.



#1 ALOS PRISM CEOS Terrain Correction problems
Hi,
The first thing to mention is that you cannot do terrain correction (a SAR process) on ALOS Prism data because Prism is optical data, not SAR. We do not yet support orthorectification of optical images. That said, I think our software should exit more gracefully with a more explanatory error message! I'll fix that!
The next error, the inability to export the DEM to a GeoTIFF has me wondering though ...Are you starting with the DEM in ASF Internal format (.img and .meta files)? I tried importing and exporting an SRTM DEM here and wasn't able to reproduce the error, so I suppose I would need a bit more information from you on how to reproduce your steps. It looks like you are getting the equivalent of a segmentation fault (the code trying to read or write to memory locations that it's not supposed to), so I'd like to look into that. If nothing else, maybe you can provide us with your DEM file that you tried to export?
Finally, the third issue is for you to produce 1B2G level data. The only difference between 1B2R and 1B2G is that the version ending in 'G' is geocoded. You should be able to geocode your Prism data with no problems by just browsing to your input file (pick the LED- file) and entering your desired settings on the Geocode and Export tabs. If you just skip the Radiometry and DEM tabs, then the Convert tool will not try to terrain correct your data and all should work fine.
Brian
#2 ALOS PRISM CEOS Terrain Correction problems
#3 ALOS PRISM CEOS Terrain Correction problems
Hi,
I'll answer your questions in the order asked:
- Orthorectification of optical images?
It has been discussed now and then, but historically the radar world has kept us pretty busy. Now that we support ALOS Avnir and Prism optical data, the topic may come up again but I can't say for sure.
- Import/Export?
I misunderstood.
- Importing Seamless Server SRTM DEMs?
Apparently my SRTM DEM had some differences than yours. When importing yours, I experienced a segmentation fault as well. We used to use separate importers (code) for various flavors of GeoTIFFs but now we are migrating everything into a single 'generic' GeoTIFF importer ...a part of the newer code (somewhat obviously) had a bug in it. It has now been fixed and will be available in the next beta or release version of the software, whichever occurs first. Thanks for pointing out the bug to us!
=> On our FTP server, you will find the imported version of your DEM. I placed it into the incoming folder under the following name:
N53E82-84g_imported.zip
When you do terrain correction and wish to use the N53E82-84g DEM, just browse to the .img version of the file instead of the TIFF and it should work fine. NOTE: Many of the Seamless Server DEMs have a lot of negative values in them ...and yours is no exception. These negative values are typically less than -900.0 and indicate "no data" for the pixels which have that value. Your file uses -32767.0 as the "no data" value. As a result, viewers may have trouble with viewing this file in a way that lets you see the topology nicely, but the ASF terrain correction already knows all about Seamless Server DEMs and automatically screens out the negative "no data" values so the terrain correction can occur properly.
=> If you have other DEMs that may be subject to the same bug, note that you can either import them with the current stable release (v3.1.16) or put them on the ASF ftp site for us to process, i.e. as I did for the N53E82-84g DEM. You can have both the stable release and the beta release installed at the same time by the way.
Regards,
Brian
#4 Re: ALOS PRISM CEOS Terrain Correction problems
Hi
I am a master student at Istanbul Tech. Univ. in Turkey. I have a thesis about generating DEM from ALOS Prism triplet mode. I have a ALOS PRISM prism triplet(ob1) dataset. Its level is 1B2R. Also I have PCI Geomatica 10.1 software.
I want to extract DEM from this dataset.Also I have GCP file(shape file format received from IKONOS at UTM projection and WGS84 datum).
PROBLEMS:
1)
I have triplet images (backward,nadir,foreword) but software doesnt support this three images at same time.I want to use all (backward,nadir,foreward) to generate one stereo image. what can ı do?
2) how many gcp is necessary for it? and extra reference DEM is necessary for it?
3)Do you have any material about generating DEM from PRISM at PCI Geomatica 10.1
note:3 images have img and vol. files and they are ceos format.
Can you help me about this matter urgently.Thanks.
Best regards.
Cihan Uysal
Geomatic Engineering
Istanbul Technical University
Istanbul/Turkey
#5 Re: ALOS PRISM CEOS Terrain Correction problems
Cihan,
uysalc wrote:I am a master student at Istanbul Tech. Univ. in Turkey. I have a thesis about generating DEM from ALOS Prism triplet mode. I have a ALOS PRISM prism triplet(ob1) dataset. Its level is 1B2R. Also I have PCI Geomatica 10.1 software.
I want to extract DEM from this dataset.Also I have GCP file(shape file format received from IKONOS at UTM projection and WGS84 datum).
PCI Orthoengine has the ability to generate a Digital Surface Model (DSM) from PRISM data. It can be a bit tricky to use, but does generate a reasonably good DSM. PCI recommends use of the L1B1 product as this is what their software is written for. A colleague of mine has been able to generate DSMs from the L1B2R products also, however.
Quote:PROBLEMS:
1)
I have triplet images (backward,nadir,foreword) but software doesnt support this three images at same time.I want to use all (backward,nadir,foreward) to generate one stereo image. what can ı do?
I am not aware of any software that deals with the three data sets of a triplet at the same time. Most packages are written to deal with stereo pairs. So I believe you need to work with your data set in pairs and combine the results later.
Quote:2) how many gcp is necessary for it? and extra reference DEM is necessary for it?
The exact number of your GCPs is not as relevant as the quality of them. You can make it work with very few.
Quote:3)Do you have any material about generating DEM from PRISM at PCI Geomatica 10.1
note:3 images have img and vol. files and they are ceos format.
I don't have any material. I did not get the time to have a closer look at the various software packages available. I have a background in photogrammetry, so I am certainly familiar with the concepts. I hope to have a chance to look into this some more.
I hope this answers your questions.
Cheers,
Rudi
#6 Re: ALOS PRISM CEOS Terrain Correction problems
Hi Cihan !
PCI recommend at least 6 GCPs to each L1b1 CCD.
Our team is studying the DEM extraction and orthorectification using also 1b2R level. But we don't have final results yet.
The GCPs must have elevation information. If you don't have it, you can extract from one reference DEM. The DEM quality has influence on the sensor modeling though.
You can download the PRISM dem extraction tutorial from PCI webpage. It will help you.
http://www.pcigeomatics.com/services/support_center/tutorials/pdf/dem_extraction_prism.pdf
regards
Aline
#7 Re: ALOS PRISM CEOS Terrain Correction problems
Hi
I am also given a new task to extract DEM form prism tripplet using Geomatica 10.0 or using LPS. I have viewed all discussions in this forum and feel it better to use geomatica. Now I have a confusion to purcahase ALOS PRISM data. It is clear I must go to level IB1 but the problem is that in alos description guide I have find two terms i.e Scene ID and Uncorrected Scene ID. I could not distinguish between these two, please tell me the difference between these two and also tell me what I will need for DEM extraction simple or uncorrected and also please inform me does geomatica 10.0 support CEOS format or I have to purchase data in different format e.g geotiff etc
thanks
Hamdam
#8 Re: ALOS PRISM CEOS Terrain Correction problems
hamdam wrote:the problem is that in alos description guide I have find two terms i.e Scene ID and Uncorrected Scene ID. I could not distinguish between these two, please tell me the difference between these two and also tell me what I will need for DEM extraction simple or uncorrected
Hi Hamdam,
what alos description guide are you referring to? Can you please provide me a link to it? Of the top of my head, I do not know what you are referring to when you say Uncorrected Scene ID.
Thanks,
Wade