FORUM FORUM
Combination View Flat View Tree View
Threads [ Previous | Next ]
RE: Co-registration of TerraSAR-X and Envisat ASAR image
toggle
Dear users,

I've been trying to co-register an Envisat ASAR IMP image (*.N1 file) with a Terrasar-X image ( TSX1_SAR__EEC_RE___SC_S_SRA_* file) using NEST 4B.

Although I have both images in the same projection WGS84 UTM 29S it is not working. During the co-registration process the TSX-Image is rotated and does not match the Envisat ASAR (which I use as the master image). What is wrong? Maybe the input into the co-registration chain? What is required as input here? ...
Interestingly, when I look at the projected input images I am told by the running text below the image that both images are in UTM 29S WGS84 ... while when I open the "co-registered" images ... it says "Satellite coordinates" ...
I also tried to re-project both images to a different projection, e.g. polarsterographic ... with the same result ... while both images look fine when I look at them separately ... when I do the co-registration the TSX-image becomes rotated again.

Perhaps somebody can help?

Thank you!

Best Regards

Stefan Kern
RE: Co-registration of TerraSAR-X and Envisat ASAR image
9/5/11 5:10 PM as a reply to Stefan Kern.
What are the residual values? Try using a much larger window size for the cross correlation in the GCP selection.
Also because these are detected images, try with and without resampling in the create stack step. (resampling here is not recommended for complex products)
RE: Co-registration of TerraSAR-X and Envisat ASAR image
9/7/11 10:22 AM as a reply to Luis Veci.
I don't think co-registration of images from different sensors will necessarily work unless the geometry of the two images is very similar (both ascending or descending, similar look angle etc.).

You can try orthorectifying the images separately with a good enough DEM - the pixels should align after that.

Marcus
RE: Co-registration of TerraSAR-X and Envisat ASAR image
9/12/11 10:35 AM as a reply to Marcus Engdahl.
Dear Marcus,

thanks a lot for this valuable comment! Now it seems more clear to me why I do not have any chance to co-register these two (ASAR and TSX) images no matter what I try to tell NEST.

However, since these scenes are from the Antarctic I wonder whether a "good enough" DEM does exist at all.

What other possibility do I have with NEST to accurately overlay these two images ... or also ALOS PALSAR and ASAR imagery. I mean, we do not live in a 1-frequency world, right?

Kind regards and thanks a lot again for your helpful comment.

Stefan Kern
RE: Co-registration of TerraSAR-X and Envisat ASAR image
9/12/11 10:41 AM as a reply to Luis Veci.
Dear Luis,

thanks a lot for your fast help. However, in the meantime, I am close to giving up because I think I have tested to co-register these two images now with quite a number of settings such as more GCPs, large or small search windows, large or small row or column interpolation factors, different accuracy thresholds and, finally, with and without sampling. It simply does not work and the images never get aligned onto each other the way they should. Either they are different in size and/or rotated against each other or they are not aligned on top of each other but at least have the same size and are not rotated. But I have no clue how to proceed further. I even tried a different projection for the ASAR ...
One idea is ... that the Enhanced Ellipsoid Correction (EEC) used for the TSX image cannot be reproduced for the ASAR image (Ellipsoid Correction possible ... but not an "enhanced" one), so maybe this is the problem?

Also, the comment of Marcus Engdahl is not too encouraging. If what he says is true with NEST then I would also not be able to overlay RADARSAT-2 and Envisat ASAR properly on top of each other .... or PALSAR and ASAR?

I would appreciate further comments.

Thank you so much!

Best regards,

Stefan Kern
RE: Co-registration of TerraSAR-X and Envisat ASAR image
9/12/11 11:23 AM as a reply to Stefan Kern.
Hi,

You are trying to coregister TSX and ASAR image of the same area, and acquired in the same orbit (ascending or descending)?

Quickly few notes and pointers you might find useful:

---
DEM over Antarctica:

you can try this one, The Radarsat-1 DEM (RAMP), http://nsidc.org/data/nsidc-0082.html . Btw: be aware that the projection of this DEM is polar polar stereographic, since neither lat/lon nor UTM makes sense there. I attached the subset of DEM to this thread.

Coregistration:
---
as long as you do not want to interferometrically combine coregistered images from different sensors/modes you should be able to coregister and do something with coregistered data. Still it depends on the application, and what exactly you want to do with the data after the coregistration.

Approach Marcus indicated is a straightforward and simplest way of getting all your images in the same reference geometry (UTM or whatever). Note that you also have to have good orbits here, and that the timing and GCP annotations would influence the overall accuracy of the end result.

Correlation optimization approach: this one could be tricky. If you try to coregister TSX vs ASAR (directly with coregistration implementation in NEST) it will less likely give optimal results. Since you are coregistering apples and oranges. Completely different resolution, scattering mechanisms, etc.... Targets that could be associated to one object in TSX image, are having completely different radiometric characteristics then in ASAR image, etc.... Perhaps you can work out something using features, like roads, but, less likely.

Moreover(!), eg 32x32 window in TSX covers X times smaller area on the ground then in ASAR. etc, etc, etc....

What you can do, and I think it should work: to either multilook TSX image to the sampling of ASAR, or oversample TSX image to resolution of ASAR and try the coregistration then. You mentioned sub-sampling, but it is not clear you tried this.

You can perhaps also geocode TSX, and ASAR, and then coregister (with correlation optimization) ASAR.geocoded to TSX.geocoded. Since TSX has much better accuracy of all the annotations.

Still, approach and success in doing the coregistration of mixed data very much depends on the application.

Hope this somehow helps.

Good luck,
Petar

EDIT: forgot the attachment.
Attachment

Attachments: RAMP_dem_subset.png (253.5k)
RE: Co-registration of TerraSAR-X and Envisat ASAR image
9/21/11 12:06 PM as a reply to Petar Marinkovic.
Dear Petar and others following up our conversation,

thanks a lot for your answer!

The reason why I want to co-register a TSX with an ASAR image

(both from decending orbits, ASAR is a stripmap IMP (swath IS2) mode image, type ASA_IMP_1P, range and azimuth spacing is 12.5 m, TSX is a level1B ScanSAR image of type EEC_RE_SC_S, range and azimuth spacing is 8.25 m)

is that I want to carry out an ice type classification with subsequent thin-ice thickness retrieval.

I really appreciate your tips and have tried an oversampling (does not work for TSX because specifying 12.5 m does not work), tried an undersampling to 12.5 m (worked) and tried a multilooking to 100 m spacing ... Trying to automatically co-register the undersampled TSX with the ASAR image worked but did not yield a satisfying result: images are shifted considerable and are even rotated a bit (2°) towards each other. Neither a different number of GCPs nor different sizes of the correlation window make a difference here. The multilooked images cannot (according to NEST) be co-registered automatically because they do not overlap (at least this is the message NEST gives).

Now, by chance I superposed both (original) images with a graticule overlay and figured out that the ASAR IMP-mode image is fine! BUT the TSX image seems not to be imported into NEST the correct way because the graticule overlay is rotated relative to what is should be (which can be seen in the preview image) by about 35 degrees ... so I suspect that the fact that the TSX image is not imported correctly creates this mess ....?

Kind regards,

Stefan Kern
RE: Co-registration of TerraSAR-X and Envisat ASAR image
9/21/11 12:48 PM as a reply to Petar Marinkovic.
Dear Petar,

I guess I've found the problem: The tie point grids (latitude, longitude, incidence angle, etc.) which are provided with the TSX data have the same size as the TSX image, yes. And these should be the latitude and longitude of each TSX-grid cell ... but obviously NEST treats these grids as if these coordinates are the coordinates of the entire shown areas, i.e. the TSX image plus the NaN values around.

In other words, if I plot the TSX-image in NEST it shows up exactly as it should as can be seen in the preview image, rotated towards the vertical with areas covered by NaNs in the 4 corners of the shown rectangle area.

If I plot the latitude, however, I get a full coverage, no NaN values; same is true with longitudes.

If I overlay the TSX-image, the latitude image and the longitude image with a graticule overlay, then this fits with the latitude and longitude image but not with the TSX-image because the latter has a different spatial coverage ... how nice ... so I have to convince NEST now that the given lat / lon values are only for the area with "good" TSX-image values ... let's see....

Kind regards,

Stefan Kern
NEST Forum Participation NEST Forum Participation

We encourage you to  join the NEST community and participate in the forum. The forum is maintained by the NEST project team who will most likely answer your technical questions related to NEST if not done so by other community members. Before posting on the Forum please have a look at the User manual, FAQ and previous posts.

Chat
Chat is temporarily unavailable.