Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1630 1367 (N07W47 L=307)   607    108        2 (1596/1334-1664/1391) -31.0 W44-W50 N05-N08
      2 1699 1359 (N04W53 L=312)   648    122        0 (1665/1327-1740/1383)   0.0 W48-W57 N02-N05
      3 1289 1301 (N13W23 L=282)   978    317        0 (1160/1217-1428/1356)   0.0 W12-W33 N11-N14
      4  342 1029 (N22E44 L=216)  3509    710        6 ( 141/ 975- 560/1132) -26.6 E62-E28 N14-N32
      5 1719  845 (S25W39 L=299)  3082    241        4 (1597/ 660-1799/1032)  33.7 W26-W50 S33-S16
      6  897  695 (S10E16 L=244)   228     69        1 ( 861/ 680- 937/ 715) -45.0 E17-E13 S11-S07
      7  995  573 (S20E14 L=245)   274     77        1 ( 959/ 552-1037/ 589) -11.3 E16-E11 S20-S18
      8  650  351 (S25E46 L=214)  2442    405        4 ( 599/ 220- 703/ 472)  57.7 E53-E40 S35-S16
      9  514  340 (S22E57 L=202)   547    175        0 ( 498/ 298- 532/ 382)  33.7 E61-E52 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