Hi everyone! I’m working with HARP in R, currently with SCIAMACHY CH4 data. I’m fairly new so it’s having a bit of a hard time fully understanding what I’m doing at HARP.
So far I have managed to obtain the CH4 data from the SCIAMACHY files (format .N1), apply a couple of quality filters and do the spatial binning in my area of interest (Argentina).
This is the code:
SCI_L3 = harp::import(file,
operations = ‘CH4_column_number_density_uncertainty > 0;
solar_zenith_angle < 80;
cloud_fraction < 0.2;
derive(CH4_column_number_density{time}[molec/m2]);
keep(latitude_bounds,longitude_bounds,
CH4_column_number_density);
bin_spatial(143,-58,0.3,45,-48,0.6)’,
‘dataset=nad_ir1_ch4’)
Then I add the CH4, latitude and longitude values into an array and convert it to raster (crs = WGS84). Everything seems fine, but when I visualize the plot in QGIS I get this:
Something tells me that the raster is rotated (maybe 90 ° AC?). I’m not quite sure I understand the arguments in bin_spatial (), does lat/lon_edge_offset represent the lower left corner of the grid?
I also don’t know why I have those gaps in the data…
I know this is a lot and not very specific, but any help you can give me would be great! Thank you