Reading and mapping OCN products - Windmap

Dear all,

I would first like to thank you all for the precious support that is provided in this forum. I have recently started working with EO, and I am doing my really first attempts with S1, aiming at mapping sea wind flow fields. Before creating this topic I’ve tried to find some similar issues, but I could not find any solutions.

I have problems processing OCN products. I would like to investigate sea wind fields, so I started taking advantage of the OCN products, which should provide already processed information about wind speed and direction. Everything seems to have a proper sense when importing files in SNAP. Wind data, already in m/s, are correctly disposed and the useful “worldwind analysis view” is working fine, as you can see in the screenshot.

So the idea was to actually use the data, importing them in QGIS. Some weird things are happening. I basically:

  1. Import the OCN product in SNAP
  2. Export: Vector Data —> owi_wind_data —> Geometry as Shape file

I import the vector layer (.shp) in a new QGIS project and that’s the result:

Geographical information of the wind flow field are basically some kind of “upside-down”. So, I go back to the SNAP tool, trying to implement some geometrical terrain correction, but unfortunately, it seems that this is not possible:

Do you have any suggestion of how actually import the OCN as vector field, in a correct manner, ensuring consistency with geographical information?

Thanks a lot in advance!

Carlo

I can confirm that exporting the wind vector as shapefile did not work for me either. As a test, I opened the table in SNAP, select all entries (Ctrl+A), then inserted them into an empty text file (here export.txt) and read it as a delimited text into QGIS.

Again, the points are not rotated with the projection of the data. The coordinates seem to be incorrect

Thanks a lot Andreas. That was a nice check. I’m kind of struggling because I realize the power of the OCN product results, but being unable to use them is kind of frustrating :slight_smile:

1 Like

definitely. I wonder if in this tutorial the same problem happens, but is ignored? https://www.youtube.com/watch?v=i7mY6Q-FCRw

I had the same thought. I don’t want to be rude or to judge other people work, but I’m feeling like this tutorial is kind of wrong. The results are showing wind vectors on the land. It should be impossible, since the measurement are clearly sea-based. Hence, I think that, as you said, it has been ignored. What a pity.

I fully agree…

1 Like

We have seen very few OCN product users on these forums - good to see that you are here, let’s try to fix the issues and define recommended processing chains for OCN product processing.

@lveci would you have any ideas?

4 Likes

@mengdahl Thanks a lot! By the way, I’ve just tried to play around my issue and apparently now things are fine. This is the process I’ve carried out:

  1. Export from SNAP the shapefile
  2. Import the .shp in QGIS
  3. From the .shp layer in QGIS, I export the correspondent .csv file
  4. I re-import the layer, now using the option “import from delimited text” as @ABraun has done before. The only difference is that I’m now importing the .cvs file that I generated in the step n.3, from QGIS, and not the copy-paste from SNAP.

The result is great now. Anyway, it’s not a smooth process, and to me it looks like a trick. I’m trying to figure out what’s the issue. Could it be something related to the reference system? That’s because, when importing from an external text, QGIS let you select the RS (and I’ve used WGS 84). Any other ideas? @ABraun Could you try with your wind field as well?

Thank you all!
Carlo

2 Likes

nice work around, I can confirm that this works

Import of SHP exported from SNAP

Export of SHP as CSV from QGIS and importing the CSV

Result

Still, I don’t get why the coordinates from the table are incorrect and then turn out to be correct after re-importing as CSV.
A direct SHP export from SNAP would be nice and I think it’s already quite nearly achieved.

2 Likes

Glad to know that it also worked with your case. Let’s see if this mystery will be somehow solved :slight_smile:

It looks like SNAP extract locations of measurements both :

  • as longitude and latitude (from the owiLon and owiLat variables provided in the OCN product)
  • as wkt geometry like POINT(-6.0791… 54.6164), but it is difficult to say how it is computed here as this location does not match the one of longitude and latitude above.

The discrepancy between the two set of position appear clearly on one of your snapshot.
image

On your first attempt, while exporting in Shapefile, SNAP is using the wkt geometry so as to record the geographic location (the wrong one), and adds the remaining data as attribute of each position. Thus the geometry is not good, while the attribute is fine.

Then in QGIS, when you export an shp layer in QGIS as CSV, it actually does not export the geometry, but only the table of attribute. However the metadata contain owiLat and owiLon that is correct (see previous point).

Then when your reopen it with QGIS, it selects the geometry definition as point coordinates associated to owiLat and owiLon, based on assumption on the names, and it is working fine.

This fully explains why the display is first incorrect then correct.
This confirms as well that there is an issue with the computation of the geometry of each point in OCN products within SNAP.

(Mistery solved)

2 Likes

thank you for this investigation and the well understandable difference.

Great! Thanks a lot :slight_smile:

I had the same problem and after using the .csv method this seems to work. However, that is quite frustrating since the “World wind analysis view” and “World Map” correctly locate the image. It seems that it comes from a wrong reading of the data in the explorer. Did you found a method to correctly apply the coordinate, not only for the vector files, but also for the wind raster?

Anymway, thank you all for your interest in the subject.

Hi @RomainJatiault. No, I did not find any other method besides the one already posted here. Thanks a lot to you for your interest. If in the meanwhile you managed to have some results with your rasters then let us know :slight_smile:

This is just to provide an update on this topic; our work based on the use of OCN products here discussed has been recently published - here the link, the paper is open access The use of Sentinel-1 OCN products for preliminary deep offshore wind energy potential estimation: A case study on Ionian sea - ScienceDirect

Cheers!
Carlo

2 Likes

Hey please can you help me find this as LAT/LON

47.552 ± 39.66 meters, -52.723 ± 39.66 meters