Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  449 1445 (N30E37 L=193)  1160    279        4 ( 359/1398- 534/1496) -28.1 E45-E29 N26-N34
      2 1042 1347 (N18W04 L=234)   272     67        1 ( 995/1327-1083/1368) -21.8 W01-W06 N16-N19
      3 1829 1300 (N09W61 L=291)   472    123        0 (1812/1255-1844/1346)  59.0 W58-W62 N06-N11
      4  298 1303 (N22E49 L=182)   491    158        2 ( 253/1264- 336/1335) -29.1 E53-E44 N19-N24
      5  301 1117 (N11E48 L=183)  1055    235        4 ( 253/1040- 347/1218) -70.0 E52-E43 N06-N17
      6 1258 1142 (N04W15 L=245)   382     90        0 (1242/1094-1275/1195) -90.0 W14-W15 N00-N06
      7 1570 1086 (S02W35 L=265)   641    135        2 (1517/1041-1634/1121) -35.5 W31-W39 S05-N00
      8  220  960 (N02E57 L=173)   938    265        6 ( 165/ 893- 272/1056)  26.6 E64-E50 N00-N07
      9  971  948 (S05E04 L=227)   239     82        0 ( 923/ 918-1008/ 982) -45.0 E06-E01 S07-S02
     10 1402  888 (S13W22 L=252)   455    102        2 (1367/ 844-1429/ 927) -53.1 W19-W23 S15-S10
     11  883  866 (S09E10 L=220)   195     34        0 ( 853/ 847- 903/ 883)  18.4 E11-E08 S10-S08
     12  210  728 (S11E63 L=167)  4666    823       16 ( 111/ 578- 296/ 871) -69.8 E77-E55 S20-S01
     13 1189  801 (S16W08 L=238)   182     17        1 (1170/ 785-1210/ 820) -26.6 W06-W09 S16-S14
     14 1264  729 (S21W13 L=243)   627    232        1 (1181/ 667-1324/ 818)  31.0 W06-W16 S24-S15
     15 1373  719 (S23W20 L=250)   694    196        2 (1309/ 632-1456/ 807)  45.0 W15-W25 S27-S17
     16  492  535 (S26E43 L=188)   874    154        2 ( 439/ 505- 551/ 569)  11.3 E47-E37 S27-S24
     17 1408  389 (S46W25 L=255)  1539    219        7 (1302/ 352-1492/ 450)  10.0 W15-W33 S49-S41
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y