Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  514 1607 (N15E51 L=262)   313     74        0 ( 499/1573- 527/1636) -45.0 E53-E48 N12-N16
      2  310 1536 (N07E67 L=246)  1718    148        2 ( 262/1472- 358/1612)  79.7 E71-E64 N01-N12
      3  925 1421 (N14E17 L=296)   322    105        0 ( 897/1379- 949/1478) -81.9 E17-E15 N10-N17
      4 1633 1105 (N17W34 L=347)   377     83        0 (1607/1071-1672/1129)  59.0 W32-W36 N14-N19
      5 1666  904 (N07W42 L=355)  2062    374        5 (1576/ 793-1714/ 983) -19.7 W33-W48 N01-N11
      6 1835  809 (N08W61 L= 14)   989    233        1 (1805/ 737-1863/ 885) -26.6 W55-W68 N04-N10
      7  541  718 (S36E23 L=290)   500    175        0 ( 481/ 685- 620/ 743)  35.0 E27-E17 S39-S32
      8 1670  443 (S17W65 L= 18)  1144    251        1 (1646/ 392-1690/ 512) -24.4 W60-W71 S18-S13
  

    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