Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1389 1756 (N57W31 L= 17)   369     59        0 (1361/1738-1418/1772)  21.8 W26-W36 N54-N58
      2  695 1620 (N37E32 L=314)  5034    681        7 ( 451/1539- 983/1765)  -7.7 E82-E08 N31-N44
      3 1055 1469 (N30E02 L=344)   375    105        0 (1026/1421-1083/1509) -56.3 E04-W00 N26-N32
      4  864 1420 (N25E14 L=332)   343    138        1 ( 846/1356- 879/1491) -77.5 E15-E13 N20-N29
      5  758 1394 (N22E21 L=325)   413    165        1 ( 738/1326- 794/1464) -53.1 E23-E17 N18-N26
      6 1497 1378 (N27W30 L= 17)   451    150        0 (1434/1334-1551/1399)  -6.3 W25-W34 N24-N28
      7 1666 1173 (N15W43 L= 29)   471    130        0 (1610/1157-1728/1199)   0.0 W38-W48 N14-N16
      8 1651  960 (N02W42 L= 28)   464    136        0 (1611/ 903-1707/ 991)  29.7 W39-W46 S01-N04
      9  275  864 (S15E52 L=294)   746    205        0 ( 219/ 823- 333/ 908)  22.6 E58-E46 S17-S11
     10  399  829 (S16E41 L=305)   257     75        1 ( 383/ 808- 424/ 865) -71.6 E42-E38 S17-S13
     11  644  570 (S30E23 L=323)   726    167        2 ( 577/ 534- 730/ 600)  26.6 E28-E16 S33-S27
     12  446  417 (S45E48 L=298)   527     55        2 ( 416/ 388- 475/ 447)  -8.1 E52-E44 S47-S42
  

    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