Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1302 1730 (N54W18 L= 83)   337     21        0 (1277/1717-1327/1744)   9.5 W15-W21 N52-N55
      2  595 1474 (N25E36 L= 29)  1575    344        2 ( 506/1394- 653/1606) -86.2 E41-E31 N19-N34
      3 1489 1407 (N30W30 L= 95)  2551    562        3 (1266/1233-1621/1530) -30.7 W12-W39 N19-N37
      4 1654 1268 (N22W43 L=108)   282     52        1 (1634/1236-1675/1293)  45.0 W40-W45 N20-N23
      5 1354 1198 (N15W19 L= 85)   437    157        1 (1322/1146-1398/1272) -54.5 W17-W22 N12-N19
      6 1537 1200 (N17W33 L= 98)   376     95        0 (1514/1143-1550/1247) -80.5 W30-W33 N13-N19
      7 1671 1189 (N17W44 L=109)   273     59        0 (1646/1172-1700/1204)  11.3 W41-W46 N16-N18
      8  438  782 (S19E38 L= 27)   735    137        2 ( 415/ 713- 473/ 843)  69.4 E40-E35 S23-S15
      9 1452  650 (S18W33 L= 98)   849    190        1 (1370/ 588-1525/ 710)  48.4 W28-W37 S22-S13
     10 1715  652 (S16W55 L=120)   538     34        0 (1695/ 622-1736/ 676)  71.6 W53-W58 S17-S14
     11 1719  568 (S21W60 L=125)  1508    157        4 (1659/ 514-1787/ 615)  49.4 W53-W69 S25-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