@@ -43,7 +43,7 @@ int spectral_predict(ard_t ard, small *mask_, int nc, int sid);
43
43
#define _SPECHOMO_N_SIM_ 10 // max. number of close clusters to be used for kNN
44
44
#define _SPECHOMO_GOOD_SAM_ 0.0698132 // good clusters are closer than this angle (4°)
45
45
#define _SPECHOMO_MEDI_SAM_ 0.2094395 // ok'ish clusters are closer than this angle (12°)
46
- #define _SPECHOMO_POOR_SAM_ 0.2617995 // clusters further away from this angle (15°) should not be uses
46
+ #define _SPECHOMO_POOR_SAM_ 0.2617995 // clusters further away from this angle (15°) should not be used
47
47
// also used to compute a weight for each cluster
48
48
#define _SPECHOMO_MIN_WEIGHT_ 1.0 - _SPECHOMO_GOOD_SAM_/_SPECHOMO_POOR_SAM_
49
49
#define _SPECHOMO_MED_WEIGHT_ 1.0 - _SPECHOMO_MEDI_SAM_/_SPECHOMO_POOR_SAM_
@@ -747,13 +747,14 @@ double pred, wpred[_SPECHOMO_N_DST_], wsum;
747
747
#ifdef FORCE_DEBUG
748
748
printf ("\n" );
749
749
printf ("on/off: %d\n" , ard .msk [p ]);
750
+ printf ("cluster number: %d\n" , s );
750
751
printf ("ard: " );
751
752
for (b = 0 ; b < _SPECHOMO_N_SRC_ ; b ++ ) printf ("%05d " , ard .dat [b_src [b ]][p ]);
752
753
printf ("\n" );
753
754
printf ("lib: " );
754
755
for (b = 0 ; b < _SPECHOMO_N_SRC_ ; b ++ ) printf ("%05d " , _SPECHOMO_CENTER_ [sid ][b ][s ]);
755
756
printf ("\n" );
756
- printf ("xx: %.2f, yy: %.2f, xy: %.2f, sam: %.5f, weight: %.5f \n" , xx , yy , xy , sam , weight [s ]);
757
+ printf ("xx: %.2f, yy: %.2f, xy: %.2f, sam: %.8f, sam_deg: %.8f, weight: %.8f \n" , xx , yy , xy , sam , sam * _R2D_CONV_ , weight [s ]);
757
758
#endif
758
759
759
760
// maximum weight of close clusters (-1 if no close cluster)
0 commit comments