Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  662 1662 (N22E41 L=239)  1507    168        4 ( 588/1612- 750/1698)  41.2 E44-E35 N17-N25
      2 1562 1493 (N35W23 L=302)  1295    224        1 (1450/1436-1656/1544)   3.2 W13-W31 N33-N37
      3  316 1506 (N06E62 L=218)   739     79        2 ( 278/1472- 343/1542)  45.0 E64-E59 N03-N08
      4  453 1440 (N04E46 L=234)   327     79        0 ( 439/1403- 466/1484) -68.2 E47-E44 N01-N06
      5  213 1297 (S08E63 L=216)  1132    224        4 ( 173/1224- 263/1390) -63.4 E68-E55 S11-S01
      6 1820 1198 (N25W49 L=329)  1064    264        1 (1792/1089-1842/1331) -86.4 W46-W52 N17-N33
      7  850 1201 (S02E14 L=266)   319     72        1 ( 809/1171- 889/1229)   0.0 E16-E11 S03-N00
      8  855  850 (S21E06 L=274)   618    179        0 ( 777/ 791- 916/ 886) -11.3 E10-W00 S22-S19
      9 1013  679 (S27W08 L=288)   201     75        0 ( 982/ 647-1036/ 708) -11.3 W05-W10 S27-S25
     10  566  550 (S47E20 L=259)  1214    262        1 ( 461/ 499- 691/ 605)  -2.5 E31-E08 S48-S43
  

    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