Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1583 1583 (N19W55 L=328)   462     89        0 (1556/1546-1606/1613)  53.1 W52-W56 N17-N21
      2  712 1469 (N40E07 L=266)   403     96        1 ( 663/1437- 768/1490)  -6.3 E11-E02 N38-N40
      3 1432 1457 (N17W37 L=310)   513     83        2 (1395/1434-1479/1490) -59.0 W34-W39 N14-N20
      4 1311 1246 (N10W22 L=295)  1903    661        5 ( 976/1136-1580/1299) -11.3 W00-W40 N02-N14
      5 1669 1207 (S03W44 L=317)   264     57        1 (1649/1181-1692/1225)  14.0 W42-W46 S03-S01
      6  775 1113 (N17E12 L=261)   362    179        1 ( 693/1084- 852/1165) -45.0 E15-E07 N13-N21
      7  291  943 (N18E50 L=223)  3015    648        5 ( 118/ 890- 494/1014) -18.4 E67-E33 N11-N24
      8 1560  752 (S25W25 L=298)  2390    458        6 (1428/ 578-1643/ 928)  32.5 W12-W34 S33-S17
      9  834  488 (S20E27 L=246)   575    129        3 ( 772/ 446- 894/ 512) -50.2 E30-E23 S24-S17
     10  479  358 (S20E59 L=214)  3845    559        5 ( 309/ 188- 606/ 506) -84.1 E68-E53 S36-S07
     11  754  424 (S22E35 L=238)   377     74        1 ( 721/ 388- 782/ 463) -90.0 E35-E32 S25-S19
  

    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