Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1183 1858 (N57E24 L=356)  1439    165        1 (1100/1814-1247/1884)   0.0 E33-E14 N54-N60
      2  518 1707 (N23E60 L=320)  2714    250        3 ( 487/1606- 557/1815) -27.3 E80-E49 N16-N32
      3 1598 1650 (N56W29 L= 49)   826     61        1 (1573/1602-1622/1686)  63.4 W25-W35 N50-N60
      4 1823 1033 (N19W54 L= 75)   710     78        2 (1795/1001-1850/1068) -26.6 W50-W57 N16-N20
      5  369  936 (S26E41 L=339)   929    278        1 ( 320/ 824- 404/1046) -50.7 E46-E36 S31-S20
      6 1884  931 (N16W65 L= 85)   514     45        2 (1867/ 905-1907/ 956)  63.4 W62-W69 N13-N17
      7  517  788 (S31E26 L=354)   305    119        0 ( 463/ 772- 571/ 812)  32.0 E30-E22 S33-S28
      8 1239  674 (S18W23 L= 43)   205     65        0 (1203/ 662-1275/ 686)  36.9 W20-W24 S19-S16
      9  582  359 (S62E15 L=  5)   519     69        0 ( 551/ 334- 622/ 387) -12.5 E19-E09 S63-S59
     10 1148  348 (S40W34 L= 54)   295     54        1 (1119/ 337-1181/ 367)  21.8 W31-W36 S41-S39
     11 1346  334 (S34W49 L= 69)   924     99        0 (1299/ 312-1403/ 361)  29.7 W43-W53 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