Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1232 1784 (N56E05 L= 35)  1826    249        3 (1115/1741-1329/1816)  31.0 E15-W05 N49-N62
      2 1152 1404 (N24W00 L= 40)   375    128        0 (1096/1378-1207/1436)  -7.1 E04-W04 N22-N24
      3 1792 1075 (N17W53 L= 93)  4127    514       10 (1683/ 943-1912/1223) -55.5 W41-W70 N10-N26
      4 1860 1105 (N21W61 L=101)   866     88        3 (1824/1074-1893/1137)  20.0 W56-W67 N19-N23
      5  319 1065 (S13E48 L=351)   401    101        0 ( 289/1015- 357/1101)  50.2 E50-E45 S15-S09
      6  913  770 (S19E02 L= 38)   426    158        0 ( 845/ 748- 962/ 806)  26.6 E05-W01 S19-S15
      7  181  751 (S35E69 L=331)   870     84        3 ( 155/ 713- 206/ 787) -51.3 E74-E64 S37-S32
      8 1366  564 (S22W33 L= 72)   709    215        2 (1299/ 523-1449/ 658)  65.6 W29-W35 S25-S14
      9 1119  481 (S33W19 L= 59)   238     51        0 (1100/ 455-1138/ 509) -26.6 W17-W21 S34-S31
     10 1512  382 (S30W56 L= 95)  2301    377        1 (1354/ 321-1657/ 456)  38.3 W45-W65 S37-S21
  

    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