2023-05-23

Problem with artifacts when ploting TROPOMI data

Hello I experienced problems of artifact when ploting TROPOIMI SO2 data using PlateCarree() projection with cartopy, loading TROPOMI data via xarray.

The problem occurs when most of the data are around -180 and 180 longitude.

The problematic figure is attached below

Image in Platcarre()

The same data are correctly plotted using Orthographic projection.

Orthographic

So I assume that the problem is with plotting itself.

Here is the code



xr_data = xr.open_dataset(file_name, group='PRODUCT', 
                              engine='netcdf4', decode_coords=True)
xr_data=xr_data.set_coords(("latitude", "longitude"))
xr_data_so2 = xr_data['sulfurdioxide_total_vertical_column'][0]
xr_data_so2=xr_data_so2 * xr_data_so2.multiplication_factor_to_convert_to_molecules_percm2
plt.figure(figsize=(14,6))
ax = plt.axes(projection = ccrs.PlateCarree())

xr_data_so2.plot.pcolormesh(ax=ax, x='longitude',y='latitude', robust=True,shading='none',cmap='jet')

I tried to print it.



No comments:

Post a Comment