Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  869 1604 (N26E26 L=223)  1894    367        4 ( 750/1487-1015/1685)  67.2 E29-E19 N15-N34
      2 1068 1586 (N30E13 L=236)  1186    191        4 (1008/1510-1117/1651) -48.4 E18-E08 N25-N34
      3 1522 1310 (N26W24 L=273)  2705    541        3 (1385/1066-1590/1481) -54.0 W11-W31 N11-N34
      4  202 1220 (S11E63 L=186)   911     91        2 ( 172/1183- 226/1273)  71.6 E66-E59 S13-S07
      5  159 1095 (S19E66 L=183)  2225    261        3 (  83/1016- 247/1190)  24.0 E84-E56 S24-S11
      6  419  956 (S22E37 L=212)   869    254        0 ( 361/ 871- 522/1009)  47.7 E41-E29 S27-S16
      7 1793  956 (N13W52 L=301)   450     89        0 (1775/ 911-1808/1005) -80.5 W51-W53 N09-N15
      8  877  835 (S18E04 L=245)   421    192        0 ( 843/ 766- 926/ 880) -81.9 E05-W00 S21-S14
      9  451  848 (S28E33 L=216)   553    130        1 ( 401/ 832- 528/ 866)  21.8 E37-E26 S29-S24
     10 1343  772 (S10W25 L=274)   395    147        0 (1279/ 730-1399/ 798)   7.1 W21-W29 S11-S08
     11  760  381 (S51W01 L=251)  1572    239        2 ( 668/ 303- 850/ 468) -15.9 E09-W12 S54-S46
     12  960  238 (S55W28 L=277)   350     11        0 ( 939/ 219- 981/ 255) -76.0 W25-W30 S57-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