Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1041 1877 (N58W48 L=328)  1800    252        2 ( 973/1793-1121/1957)   6.6 W31-W83 N52-N59
      2  644 1647 (N55E08 L=272)   533     96        1 ( 605/1615- 687/1681)  11.3 E13-E03 N53-N57
      3 1619 1521 (N15W55 L=335)   366     33        0 (1596/1498-1640/1541) -90.0 W53-W56 N13-N16
      4  613 1430 (N40E16 L=264)   237     51        0 ( 581/1417- 642/1443) -11.3 E18-E13 N39-N41
      5 1806  984 (S19W50 L=330)   492    164        0 (1780/ 906-1840/1047)  37.9 W46-W55 S23-S16
      6 1858  963 (S23W56 L=336)   536    109        0 (1832/ 914-1888/1014)  26.6 W52-W61 S24-S20
      7 1437  700 (S24W16 L=296)   399    162        0 (1367/ 658-1508/ 734)   0.0 W10-W21 S24-S23
      8 1249  495 (S31E02 L=278)   696    257        1 (1174/ 465-1332/ 576) -26.6 E06-W04 S35-S24
      9 1614  457 (S49W28 L=307)   920    131        2 (1580/ 421-1660/ 510)   6.3 W23-W33 S52-S43
     10  911  413 (S27E25 L=254)   612    114        2 ( 885/ 360- 944/ 466)  40.6 E28-E21 S30-S23
  

    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