Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  874 1746 (N53E17 L=158)   536     73        2 ( 842/1722- 915/1765)  26.6 E20-E13 N50-N54
      2  361 1535 (N33E59 L=116)  2190    210        5 ( 299/1471- 436/1617) -11.9 E87-E48 N29-N37
      3 1426 1526 (N36W30 L=206)   480     93        0 (1401/1482-1447/1573)  54.5 W27-W32 N32-N39
      4 1352 1553 (N38W25 L=200)   434     99        0 (1313/1535-1397/1572)   8.1 W21-W28 N36-N39
      5  628 1458 (N29E30 L=145)  1426    314        2 ( 487/1425- 768/1489)   0.0 E41-E18 N27-N31
      6 1460 1463 (N31W32 L=207)   274     43        0 (1436/1443-1480/1484)  26.6 W29-W33 N29-N32
      7 1475 1335 (N23W31 L=206)   375    106        0 (1459/1290-1490/1385)  74.1 W28-W31 N19-N26
      8  466 1132 (N08E37 L=138)   481    147        0 ( 440/1063- 496/1196) -54.5 E39-E34 N04-N11
      9 1800 1160 (N11W57 L=232)   647    104        0 (1756/1116-1832/1192) -35.5 W52-W60 N08-N13
     10  813 1041 (N04E13 L=162)   268     88        0 ( 796/ 997- 840/1074)  63.4 E14-E11 N01-N05
     11 1851  976 (N00W62 L=237)   584     97        0 (1836/ 933-1865/1028)  68.2 W59-W63 S02-N03
     12  433  749 (S16E40 L=135)   368     35        0 ( 415/ 678- 464/ 810) -90.0 E41-E38 S20-S11
     13 1771  751 (S14W56 L=231)  1133    119        2 (1729/ 704-1817/ 803) -20.0 W50-W62 S16-S10
     14  388  657 (S22E46 L=129)   289      0        0 ( 351/ 643- 438/ 677)   6.3 E49-E40 S22-S20
     15 1089  563 (S26W05 L=180)  1131    241        3 (1019/ 514-1171/ 649) -41.6 W00-W11 S29-S20
     16  884  564 (S26E09 L=167)   422    104        0 ( 831/ 539- 930/ 599)  23.2 E12-E05 S28-S23
     17 1376  214 (S56W46 L=221)   600     65        0 (1351/ 187-1392/ 247) -49.4 W40-W51 S59-S52
  

    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