Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1361 1464 (N32W22 L=207)  1924    497        2 (1222/1317-1542/1566) -32.5 W11-W35 N23-N39
      2  677 1394 (N23E26 L=159)   677    193        1 ( 639/1318- 736/1475) -45.0 E31-E20 N18-N28
      3 1793 1286 (N22W59 L=245)   514    133        3 (1759/1219-1820/1334)  37.9 W54-W65 N17-N25
      4 1661 1260 (N20W44 L=230)   751    134        3 (1606/1237-1710/1306)  -5.7 W39-W49 N18-N23
      5 1724  823 (S07W51 L=236)  1256    195        2 (1650/ 753-1794/ 873)  32.5 W44-W57 S11-S03
      6 1310  713 (S16W20 L=206)   335    118        3 (1276/ 674-1336/ 765)  59.0 W18-W22 S18-S12
      7 1808  681 (S16W65 L=251)   561     87        0 (1786/ 641-1839/ 716)  59.0 W63-W68 S18-S13
      8 1475  630 (S20W34 L=219)   497    111        2 (1435/ 591-1511/ 680)  56.3 W31-W35 S22-S16
      9  556  591 (S29E31 L=154)  1272    343        3 ( 424/ 531- 699/ 676)  21.0 E45-E19 S33-S21
     10 1619  635 (S19W46 L=231)   350     57        1 (1589/ 609-1643/ 659)  45.0 W43-W47 S21-S17
  

    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