Changeset 11237

Show
Ignore:
Timestamp:
04/10/07 15:21:48 (2 years ago)
Author:
dnadeau
Message:

Add threshold of 0.1 degree in latitude to allow projection
for Gaussian gridded data.

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • trunk/gdal/frmts/netcdf/netcdfdataset.cpp

    r11175 r11237  
    10331033                                      pdfXCoord[xdim-1]) * 1000); 
    10341034         
    1035         if( ( abs( nSpacingBegin ) == abs( nSpacingLast )) && 
    1036             ( abs( nSpacingBegin ) == abs( nSpacingMiddle )) && 
    1037             ( abs( nSpacingMiddle ) == abs( nSpacingLast )) ) { 
     1035        if( ( abs( nSpacingBegin ) ==  abs( nSpacingLast )     ) && 
     1036            ( abs( nSpacingBegin ) ==  abs( nSpacingMiddle )   ) && 
     1037            ( abs( nSpacingMiddle ) == abs( nSpacingLast )     ) ) { 
    10381038 
    10391039/* -------------------------------------------------------------------- */ 
    …… 
    10501050                                                pdfYCoord[ydim-1]) *  
    10511051                                               1000); 
    1052              
    1053  
    1054             if( ( abs( nSpacingBegin ) == abs( nSpacingLast )) && 
    1055                 ( abs( nSpacingBegin ) == abs( nSpacingMiddle )) && 
    1056                 ( abs( nSpacingMiddle ) == abs( nSpacingLast )) ) { 
     1052 
     1053                     
     1054/* -------------------------------------------------------------------- */ 
     1055/*   For Latitude  we allow an error of 0.1 degrees for gaussion        */ 
     1056/*   gridding                                                           */ 
     1057/* -------------------------------------------------------------------- */ 
     1058 
     1059            if((( abs( abs(nSpacingBegin) - abs(nSpacingLast) ) )   < 100 ) && 
     1060               (( abs( abs(nSpacingBegin) -  abs(nSpacingMiddle) ) ) < 100 ) && 
     1061               (( abs( abs(nSpacingMiddle) - abs(nSpacingLast) ) )   < 100) ) { 
     1062 
     1063                if( ( abs( nSpacingBegin )  !=  abs( nSpacingLast )     )  || 
     1064                    ( abs( nSpacingBegin )  !=  abs( nSpacingMiddle )   ) || 
     1065                    ( abs( nSpacingMiddle ) !=  abs( nSpacingLast )     ) ) { 
     1066                     
     1067                    CPLError(CE_Warning, 1,"Latitude grid not spaced evenly.\nSeting projection for grid spacing is within 0.1 degrees threshold.\n"); 
     1068 
     1069                } 
    10571070/* -------------------------------------------------------------------- */ 
    10581071/*      We have gridded data s we can set the Gereferencing info.       */