Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  619 1411 (N13E34 L=102)   397    115        1 ( 591/1373- 659/1461) -29.7 E37-E30 N10-N15
      2  506 1211 (S01E36 L=100)   562    184        1 ( 463/1124- 563/1272) -31.0 E40-E30 S05-N01
      3 1757 1148 (N21W47 L=183)   906    151        3 (1699/1089-1792/1203) -26.6 W41-W51 N18-N24
      4  156 1081 (S15E66 L= 70)  1607    255        1 ( 142/ 937- 176/1182) -82.4 E70-E63 S24-S09
      5 1227 1139 (N09W09 L=146)   307    105        0 (1166/1109-1275/1165)   0.0 W05-W12 N07-N09
      6  193 1055 (S16E60 L= 76)   564     78        3 ( 175/1006- 210/1102) -63.4 E62-E57 S19-S13
      7  488  918 (S19E32 L=104)   539    221        2 ( 403/ 859- 558/ 989)  47.7 E37-E27 S24-S13
      8  796  805 (S19E09 L=127)   197     99        0 ( 765/ 752- 825/ 840)  74.1 E10-E07 S23-S16
      9  605  818 (S23E22 L=114)   237     51        0 ( 571/ 802- 632/ 836)   0.0 E24-E20 S23-S21
     10  462  763 (S29E32 L=104)   553    162        0 ( 409/ 708- 505/ 817)  52.1 E36-E29 S33-S24
     11 1352  742 (S11W26 L=162)   250    109        1 (1296/ 725-1400/ 762)  33.7 W22-W28 S13-S09
     12 1135  413 (S37W24 L=160)   663    152        2 (1047/ 397-1212/ 427)  24.4 W17-W29 S39-S34
     13 1557  377 (S29W57 L=193)   831    107        0 (1507/ 353-1595/ 412)  59.0 W53-W59 S30-S25
  

    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