Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  858 1863 (N44E47 L=258)   653     57        1 ( 825/1833- 890/1886)  45.0 E50-E44 N41-N46
      2  740 1746 (N32E44 L=261)  1271    157        2 ( 692/1706- 810/1785)  -4.1 E50-E36 N28-N33
      3  420 1599 (N15E59 L=246)  2593    241        9 ( 370/1508- 484/1707) -66.8 E67-E53 N08-N22
      4  770 1568 (N21E31 L=274)   453     80        0 ( 751/1529- 791/1608) -45.0 E33-E27 N18-N23
      5 1589 1309 (N27W28 L=333)  1202    168        0 (1538/1174-1631/1436) -54.5 W22-W33 N20-N34
      6 1821 1040 (N19W53 L=358)  3707    582        6 (1727/ 821-1898/1280) -47.6 W40-W65 N07-N31
      7  329 1228 (S09E50 L=255)   329    109        1 ( 309/1179- 366/1278)  80.5 E51-E48 S11-S05
      8 1672 1172 (N22W37 L=342)   415      0        0 (1628/1141-1694/1228)  66.8 W33-W38 N18-N25
      9 1297 1085 (N06W14 L=319)   475    152        2 (1258/1020-1337/1159)  90.0 W12-W14 N01-N10
     10  193 1070 (S21E61 L=244)   367     67        2 ( 174/1044- 224/1103) -26.6 E63-E56 S21-S18
     11 1766  984 (N14W49 L=353)   319    103        1 (1740/ 942-1785/1011)  90.0 W46-W50 N11-N15
     12  672  916 (S19E18 L=287)   238      4        0 ( 651/ 896- 692/ 935)  45.0 E19-E16 S20-S17
     13 1744  779 (N02W53 L=358)   935    157        1 (1680/ 749-1807/ 824)   8.1 W45-W59 N00-N04
     14  723  710 (S30E10 L=295)   274    135        0 ( 663/ 667- 773/ 730)  40.6 E13-E06 S34-S27
     15  299  603 (S49E51 L=254)  2283    225        4 ( 211/ 542- 397/ 669)  -9.1 E67-E37 S52-S45
     16 1338  617 (S18W30 L=335)   274     88        0 (1301/ 571-1369/ 651)  90.0 W29-W31 S21-S15
     17 1465  616 (S15W38 L=343)   525    145        0 (1400/ 603-1543/ 631)  24.0 W34-W43 S16-S12
     18 1224  273 (S42W45 L=349)  1608    205        2 (1155/ 232-1330/ 322)  10.6 W35-W54 S45-S38
  

    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