Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  841 1848 (N63W22 L= 39)   605     51        2 ( 811/1819- 868/1876)  26.6 W17-W27 N60-N65
      2  710 1790 (N64W02 L= 19)  1434    194        2 ( 632/1737- 790/1841)   9.5 E09-W15 N59-N66
      3  665 1476 (N40E12 L=  6)  1129    303        2 ( 533/1434- 807/1515) -26.6 E22-E02 N33-N44
      4 1631 1404 (N09W48 L= 65)   947    192        1 (1580/1330-1675/1491)  78.7 W45-W50 N04-N14
      5 1430 1402 (N15W34 L= 51)   374     81        1 (1388/1383-1472/1425)  -9.5 W30-W36 N13-N15
      6  849 1336 (N26E02 L= 15)   491    117        0 ( 802/1289- 897/1373)   9.5 E05-W01 N23-N27
      7  187 1336 (N42E64 L=313)   501     53        0 ( 158/1323- 220/1351) -11.3 E72-E57 N40-N43
      8 1566 1296 (N05W39 L= 56)   232     78        0 (1528/1283-1606/1316) -36.9 W37-W41 N03-N07
      9 1832 1266 (S04W61 L= 78)   367     11        0 (1814/1244-1849/1285) -90.0 W59-W62 S04-S02
     10 1573 1166 (S02W36 L= 53)   355     85        0 (1546/1124-1596/1204)  26.6 W32-W38 S03-N00
     11  154  819 (N11E69 L=309)  1506    153        3 ( 103/ 769- 215/ 870) -26.6 E77-E61 N07-N14
     12  762  558 (S16E27 L=350)   801    208        2 ( 709/ 460- 810/ 663)  37.6 E34-E21 S20-S10
     13 1725  563 (S44W41 L= 58)   484     87        1 (1686/ 547-1768/ 587) -10.3 W36-W47 S46-S40
     14 1531  483 (S42W21 L= 38)   883    170        1 (1444/ 465-1609/ 512) -21.8 W13-W28 S45-S39
     15  315  462 (S13E69 L=309)   721     73        1 ( 291/ 427- 331/ 501)  29.7 E73-E63 S14-S10
     16  455  457 (S15E54 L=324)   438     75        0 ( 414/ 444- 495/ 471) -15.9 E56-E49 S16-S13
     17 1309  387 (S42W01 L= 18)   505     97        1 (1269/ 355-1360/ 415)  14.0 E03-W05 S43-S40
  

    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