Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1182 1858 (N57E25 L=356)  1458    177        0 (1093/1814-1245/1883)   3.0 E34-E15 N54-N60
      2  515 1704 (N23E60 L=321)  2503    255        3 ( 488/1609- 555/1814) -30.7 E78-E49 N16-N32
      3 1594 1653 (N56W28 L= 49)   841     89        1 (1568/1604-1619/1694)  65.6 W24-W34 N50-N61
      4 1561 1306 (N27W26 L= 47)   291    103        0 (1515/1269-1601/1335)  -8.1 W22-W29 N25-N27
      5 1815 1037 (N19W53 L= 74)   814     70        3 (1789/1002-1847/1073) -31.0 W49-W57 N16-N20
      6  366  931 (S26E41 L=340)   918    310        0 ( 316/ 811- 397/1048) -55.3 E47-E37 S32-S19
      7 1883  949 (N17W64 L= 85)   585     94        1 (1865/ 906-1898/ 998)  80.5 W61-W66 N13-N19
      8  511  790 (S31E27 L=354)   306    119        0 ( 457/ 773- 566/ 814)  32.0 E31-E23 S33-S28
      9  657  730 (S31E15 L=  6)   281     76        0 ( 624/ 715- 696/ 748)  31.0 E17-E12 S32-S28
     10  579  360 (S62E16 L=  5)   451     65        1 ( 549/ 334- 616/ 387)  -6.3 E20-E10 S63-S59
     11 1144  350 (S40W33 L= 54)   299     60        2 (1114/ 338-1182/ 371)  21.8 W31-W36 S41-S39
     12 1344  337 (S34W48 L= 69)   928    106        1 (1297/ 316-1400/ 364)  33.7 W43-W52 S36-S31
  

    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