Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1149 1697 (N37E14 L=298)   476     93        1 (1099/1670-1199/1715)  14.0 E18-E10 N35-N38
      2  440 1545 (N09E53 L=260)   482     64        0 ( 415/1512- 468/1577) -23.2 E56-E49 N07-N10
      3 1630 1492 (N40W27 L=340)   843    186        2 (1596/1402-1659/1574) -67.4 W23-W30 N34-N46
      4 1582 1288 (N25W26 L=339)   402    165        0 (1517/1237-1614/1322)   0.0 W21-W28 N21-N26
      5 1798 1133 (N23W47 L=360)   569    125        0 (1756/1071-1835/1171) -18.4 W42-W52 N20-N24
      6 1734  975 (N12W45 L=357)   381     46        1 (1703/ 957-1759/1001)  26.6 W42-W46 N10-N13
      7  772  915 (S19E11 L=301)   212     57        0 ( 742/ 892- 802/ 933)   0.0 E13-E09 S19-S18
      8 1653  882 (N04W41 L=354)   252     67        0 (1613/ 871-1688/ 897)  31.0 W38-W43 N02-N05
      9 1249  709 (S17W22 L=334)   506    157        2 (1180/ 678-1326/ 732)  40.6 W18-W26 S20-S14
     10  469  582 (S48E28 L=284)   362     49        2 ( 444/ 552- 490/ 610) -56.3 E30-E25 S49-S45
     11 1508  389 (S25W56 L=  8)   519     94        0 (1470/ 360-1547/ 413)   0.0 W51-W62 S26-S24
     12  931  248 (S56W29 L=342)   449     80        1 ( 906/ 216- 945/ 288) -45.0 W24-W33 S59-S53
  

    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