Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  988 1399 (N27W09 L=119)   408    153        1 ( 912/1369-1068/1439) -10.3 W03-W14 N26-N28
      2  727 1408 (N35E08 L=103)   568    150        0 ( 659/1385- 801/1425) -21.8 E12-E02 N32-N37
      3 1394 1385 (N15W32 L=143)   529    187        0 (1306/1371-1483/1403) -24.4 W26-W37 N12-N17
      4 1814 1330 (S01W62 L=173)  1037    183        1 (1761/1300-1863/1360) -20.0 W57-W68 S04-N01
      5 1216 1276 (N14W18 L=128)   588    121        3 (1169/1244-1277/1308)   0.0 W14-W21 N11-N16
      6 1749 1113 (S10W47 L=158)  2009    440        7 (1709/ 923-1822/1277)  57.7 W42-W59 S21-S02
      7  197  950 (N20E58 L= 52)   420     83        1 ( 175/ 925- 220/ 983) -53.1 E60-E55 N18-N22
      8 1489  867 (S16W23 L=133)   632    201        1 (1406/ 819-1554/ 913) -52.1 W18-W25 S20-S11
      9  712  697 (S05E26 L= 84)   254     56        1 ( 678/ 682- 745/ 711) -26.6 E28-E24 S06-S03
     10  625  362 (S23E46 L= 65)   906    229        1 ( 600/ 290- 653/ 438)  55.0 E50-E41 S28-S18
     11 1296  309 (S48E05 L=106)  3075    411       10 (1149/ 239-1402/ 384)   4.8 E19-W05 S52-S42
  

    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