How to extract a diagonal band of data from a square matrix, to reduce array size?

1 view (last 30 days)
Robert on 18 Jul 2017
Edited: David Goodmanson on 18 Jul 2017
Hi,
I am looking to extract data for the region Latitudes (North=29 and South= 22.5) and Longitudes (West= -179 and East= -160.5), which is the area where the Northwest Hawaiian Islands are located. However, as the data comes at 1km pixel resolution, I end with a Matrix=(650,1500,time), where time is 6000, so having trouble in Matlab to produce this matrix.
As the islands are distributed in a diagonal position (within this latitude and longitude box), I am wondering how could I extract just the diagonal portion, or data close to the islands in this case, and skip extracting the other data (southwest and northwest corners) to reduce the size of my final matrix. All this to be able later to analyze this data in Matlab.
Suggestions are appreciated! Thanks!
David Goodmanson on 18 Jul 2017
Hi Robert, lots of possibilities here. First, if your data is in floating point, would there be enough resolution if you converted it with a linear transformation to uint16 (0-65535), for example? That would be worth a factor of 4. If you are only interested in the islands themselves, it would not be hard to extract equal-sized blocks and concatenate them to make a 4x3 array of blocks or similar, like stamps. Or stack them in the third dimension of a 3d array. Time would be the next dimension after that.

Walter Roberson on 18 Jul 2017