Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  995 1846 (N59W38 L=290)   660     92        3 ( 964/1806-1028/1879)  15.9 W31-W46 N56-N61
      2 1391 1730 (N36W53 L=305)   651     67        2 (1356/1703-1424/1755)   0.0 W49-W58 N34-N37
      3 1652 1377 (N09W50 L=301)   406     72        0 (1611/1364-1693/1390) -26.6 W46-W53 N07-N10
      4 1804 1227 (S04W57 L=309)  1046    128        2 (1778/1158-1830/1287)  32.0 W52-W61 S07-N00
      5  874 1241 (N23E03 L=248)   267     77        0 ( 838/1213- 904/1265) -68.2 E05-E01 N20-N25
      6  842 1155 (N19E07 L=244)   658    272        0 ( 734/1107- 957/1186) -12.1 E14-W00 N15-N20
      7  878  988 (N09E09 L=243)   269    158        0 ( 848/ 917- 915/1056)  40.6 E12-E05 N05-N11
      8 1412 1008 (S04W21 L=273)   334     69        0 (1382/ 973-1434/1033)  76.0 W19-W23 S05-S01
      9  272  914 (N16E54 L=198)   427     82        1 ( 247/ 869- 293/ 945) -59.0 E55-E51 N13-N18
     10  202  872 (N15E62 L=189)   414     71        0 ( 190/ 837- 215/ 908)  45.0 E64-E59 N12-N16
     11 1747  805 (S26W42 L=294)   681    109        1 (1704/ 764-1789/ 840)   7.1 W37-W46 S27-S23
     12  160  760 (N08E73 L=179)  1069    103        2 ( 137/ 706- 185/ 799) -63.4 E77-E68 N04-N10
     13  985  478 (S24E17 L=235)   693    201        1 ( 917/ 419-1043/ 564) -74.7 E18-E15 S28-S16
     14  319  509 (S09E65 L=187)   915    110        3 ( 292/ 457- 345/ 557)  54.5 E68-E61 S12-S05
     15  479  435 (S16E53 L=198)  3749    286        7 ( 384/ 302- 558/ 550)  34.8 E66-E43 S25-S09
     16  643  431 (S19E41 L=211)  1271    196        8 ( 584/ 351- 695/ 507)  39.3 E48-E34 S24-S14
     17 1262  342 (S42E04 L=247)  2979    593       10 (1042/ 207-1505/ 504)  10.3 E26-W18 S47-S36
  

    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