dvalters / NIMROD-toolbox

Python tools for manipulating Met Office NIMROD data, creating timeseries from it, converting to ascii etc.
4 stars 3 forks source link

Upscaling issue with Mexico data (Freya) #4

Open dvalters opened 6 years ago

dvalters commented 6 years ago

Upscaling cropped radar dem does not produce the correct hydroindex. (Cells are misaligned)

Output at this stage is correct:

pastedgraphic-1

But then the following stage where it prints out an image of the upscaled hydroindex is wrong - the extents are correct but the zoom function is not scaling the image correctly (There should be a very thin band of hydroindex cells at the top of the hydroindex file, where they are just within the bounds of the radar in the upper image, but instead they are distorted and not proportional to the region where the top cells overlay the DEM.

dvalters commented 6 years ago

Thanks Freya for reporting this:

I’ve also commented out the convert_OSGB36_to_UTM30() function that was still in the timeseries script as it isn’t needed and I was worried it might project incorrectly. The inputs/outputs are attached. The cropped radar test file is printing the correct precipitation values (you can see the colours are the same), but its location is starting from the bottom-left corner of the DEM instead of across the centre of the DEM.

screen shot 2018-06-12 at 6 57 53 pm

So I’m pretty sure this is still similar to the results we were getting last week unfortunately. The printout of the hydroindex also shows an additional problem that will be the next thing to address after the timeseries and cropped radar issue. The correct number of 12 hydroindex cells is written (we had 9 last week as I had originally cropped the DEM to make it a clean rectangle, but I haven’t changed the extents of this one). But as you can see the cells are not the same size as the original radar, nor are they aligned to the original radar or the cropped radar. I think this is also similar to the results of last week.

screen shot 2018-06-12 at 6 59 03 pm