pull/13/head
parent
f0b0156af3
commit
d0651791ea
|
@ -4188,6 +4188,7 @@ void CreateSARIntensityByLookTable(QString IntensityRasterPath,
|
|||
long blockYSize = Memory1GB / looktableds.width / 8 * 2;
|
||||
|
||||
Eigen::MatrixXd SARData = SARIntensity.getData(0, 0, SARIntensity.height, SARIntensity.width, 1);
|
||||
SARData = SARData.array() * 0;
|
||||
// 分块处理
|
||||
for (int yOff = 0; yOff < looktableds.height; yOff += blockYSize)
|
||||
{
|
||||
|
@ -4199,12 +4200,11 @@ void CreateSARIntensityByLookTable(QString IntensityRasterPath,
|
|||
|
||||
for (long i = 0; i < rowData.rows(); i++) {
|
||||
for (long j = 0; j < rowData.cols(); j++) {
|
||||
const long r =round( rowData(i,j))-min_rid;
|
||||
const long c = round(colData(i, j))-min_cid;
|
||||
long r =round( rowData(i,j))-min_rid;
|
||||
long c = round(colData(i, j))-min_cid;
|
||||
|
||||
if (r >= 0 && r < SARIntensity.height && c >= 0 && c <= SARIntensity.width) {
|
||||
if (r >= 0 && r < SARIntensity.height && c >= 0 && c < SARIntensity.width) {
|
||||
SARData(r, c) = SARData(r, c) + geoData(i, j);
|
||||
|
||||
}
|
||||
}
|
||||
}
|
||||
|
|
Loading…
Reference in New Issue