Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1509 1641 (N25W52 L=210)  4317    748        7 (1268/1578-1717/1684) -20.7 W33-W71 N16-N32
      2 1592 1490 (N14W50 L=208)   617    147        0 (1534/1418-1659/1536) -68.2 W46-W51 N08-N18
      3  280 1357 (N43E47 L=111)   429     71        0 ( 264/1325- 298/1383) -59.0 E49-E45 N40-N45
      4 1838  858 (S28W51 L=209)  3838    598        5 (1697/ 665-1902/1113)  38.9 W36-W62 S36-S14
      5  204  975 (N20E57 L=101)   785    160        0 ( 179/ 905- 233/1052)  48.8 E60-E52 N16-N24
      6  863  972 (N06E10 L=148)   469    156        0 ( 787/ 929- 945/ 996)  -5.7 E15-E05 N03-N07
      7 1705  824 (S26W38 L=197)   624    167        0 (1660/ 748-1754/ 909)  35.0 W33-W43 S29-S21
      8 1189  673 (S19W00 L=159)   315    111        0 (1151/ 622-1208/ 731)  45.0 E03-W02 S21-S16
      9  767  659 (S10E24 L=134)   300    124        1 ( 755/ 592- 783/ 709)  66.8 E26-E21 S13-S06
     10 1319  519 (S33W04 L=163)   370    122        0 (1264/ 472-1375/ 577)  24.0 W00-W09 S35-S30
     11  985  446 (S28E19 L=139)   693    241        1 ( 928/ 368-1021/ 552)  34.7 E25-E12 S32-S22
     12  497  412 (S18E53 L=105)  1503    397        5 ( 444/ 286- 542/ 526)  53.1 E59-E45 S27-S11
  

    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