Atmospheric Toolbox

About the HARP category

Discussions about the HARP software: support, examples, development ideas, etc.

Hello
I want to export a NetCDF file to raster to use it on Qgis or Arcmap…

File: (S5P_NRTI_L2__NO2____20191217T181608_20191217T182108_11284_01_010302_20191217T190947.nc )

Id already installed Anaconda, Visan and HARP. (Windows64)

When i open the environment with pyton and put the “harpconvert -a bin_spatial” command but always have syntax error.

Help please!!

Note that bin_spatial is a function that requires arguments.
You will need to give the information on the target grid you want.

There is actually a whole sequence of operations that you might want to perform in order to efficiently convert the product to a file that can be used in qgis. This includes applying the quality filter, removing unneeded variables, etc.

An example convert for the file that you provided (i.e. the parameters to bin_spatial are tailored for this specific O3 file) is:

harpconvert -a 'tropospheric_NO2_column_number_density_validity>75;keep(latitude_bounds,longitude_bounds,tropospheric_NO2_column_number_density);bin_spatial(2400,-37,0.01,3300,-83,0.01);derive(tropospheric_NO2_column_number_density [Pmolec/cm2]);squash(time, (latitude_bounds,longitude_bounds));derive(latitude {latitude});derive(longitude {longitude});exclude(latitude_bounds,longitude_bounds,count,weight)' S5P_NRTI_L2__NO2____20191217T181608_20191217T182108_11284_01_010302_20191217T190947.nc o3grid.nc

If you use a different product then modify parameters to the bin_spatial function accordingly.

See also the other forum threads such as Map Sentinel 5p CH4, SO2 and CO Data (Visan) for more details.

Hello svniemeijer
Verry thank you for the explanation.
Now i can export S5P data with HARP.

I have another question: I want to export the same S5P NO2 Data with the WPLOT Visan configuration to handle on Qgis.

This configuration, same size:

Can u help me with this bin_spatial function please
Thank u for your explanation again.

Regards

I am not entirely sure what you mean. Do you mean a plot without the filtering and unit conversion?

In that case use:

harpconvert -a 'keep(latitude_bounds,longitude_bounds,tropospheric_NO2_column_number_density);bin_spatial(2400,-37,0.01,3300,-83,0.01);squash(time, (latitude_bounds,longitude_bounds));derive(latitude {latitude});derive(longitude {longitude});exclude(latitude_bounds,longitude_bounds,count,weight)' S5P_NRTI_L2__NO2____20191217T181608_20191217T182108_11284_01_010302_20191217T190947.nc o3grid.nc

Thank you a lot !

Regards

Hi I need obtain puntual information about sentinel 5p, and I’m used a similar code but I don’t be what are this numbers (2400,-37,0.01,3300,-83,0.01).
You can help my.

The meaning of the values are explained in the HARP documentation for the bin_spatial operation.
The parameters are: lat_edge_length, lat_edge_offset, lat_edge_step, lon_edge_length, lon_edge_offset, and lon_edge_step.