Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  921 1902 (N48E51 L=220)   787     74        1 ( 891/1864- 957/1921) -49.4 E55-E45 N43-N50
      2  809 1845 (N39E50 L=221)  1190    151        2 ( 774/1780- 854/1895) -47.5 E56-E43 N32-N44
      3 1107 1630 (N32E15 L=256)   655    121        5 (1075/1566-1135/1689) -68.2 E16-E11 N27-N37
      4 1353 1550 (N35W04 L=275)  1046    299        2 (1271/1421-1412/1655) -38.2 E04-W11 N27-N41
      5 1662 1579 (N51W33 L=303)   479     59        1 (1644/1540-1682/1606)  63.4 W29-W36 N46-N54
      6 1492 1172 (N16W23 L=293)   839    211        3 (1441/1096-1534/1280) -45.0 W18-W27 N12-N21
      7  496 1168 (S12E35 L=235)   364    109        1 ( 469/1122- 527/1209) -31.0 E37-E32 S13-S09
      8 1471  991 (N05W26 L=297)   462    157        2 (1417/ 931-1519/1051)  71.6 W23-W27 N00-N09
      9  758 1004 (S14E14 L=256)   198    118        0 ( 732/ 960- 799/1040)  80.5 E15-E12 S17-S11
     10 1070  642 (S26W15 L=285)   206    136        1 (1017/ 617-1132/ 662)  50.2 W11-W17 S28-S22
     11  992  407 (S42W21 L=292)   602    216        0 ( 890/ 389-1095/ 432)  18.4 W13-W28 S45-S39
     12 1306  402 (S31W41 L=312)   536     99        1 (1252/ 382-1354/ 431)  14.0 W36-W45 S33-S29
  

    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