Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1353 1649 (N45W04 L=169)   359     85        0 (1310/1633-1393/1666)  15.9 W00-W07 N42-N46
      2  480 1610 (N17E55 L=110)   845    101        5 ( 449/1577- 529/1646)  -6.3 E59-E50 N15-N19
      3  455 1544 (N13E53 L=112)   557     83        4 ( 412/1516- 483/1573)  36.9 E56-E50 N10-N14
      4 1418 1551 (N39W11 L=177)   405     81        0 (1378/1527-1454/1569)   0.0 W08-W14 N38-N40
      5 1652 1301 (N30W34 L=200)  1887    379        2 (1512/1189-1782/1411) -12.5 W21-W48 N26-N32
      6 1081 1075 (N00W02 L=167)   446    136        0 (1045/1019-1114/1133) -26.6 E01-W05 S01-N03
      7  189 1105 (S18E63 L=102)   413     66        0 ( 168/1073- 217/1130)  38.7 E66-E60 S20-S16
      8 1116  953 (S05W07 L=172)   583     96        2 (1074/ 913-1152/ 995)  80.5 W04-W08 S08-S02
      9 1026  834 (S14W05 L=170)   341     64        0 ( 993/ 814-1062/ 853)  53.1 W03-W06 S16-S12
     10 1444  556 (S19W40 L=206)   417     66        0 (1400/ 534-1479/ 571)  21.8 W36-W43 S19-S17
  

    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