Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  497 1666 (N58E46 L= 28)   658     83        0 ( 465/1638- 539/1689)  11.3 E53-E39 N56-N60
      2 1704 1455 (N19W60 L=134)   519    114        0 (1651/1420-1742/1484) -39.8 W55-W62 N16-N21
      3  729 1435 (N36E14 L= 60)   462    146        1 ( 660/1419- 802/1446)  -5.2 E19-E08 N34-N38
      4  759 1181 (N20E14 L= 60)   355    110        0 ( 716/1156- 802/1199) -26.6 E17-E11 N17-N21
      5  378 1069 (N17E44 L= 31)   562     88        1 ( 345/1023- 409/1108) -45.0 E46-E40 N14-N19
      6 1339  589 (S26W14 L= 88)  1878    651        3 (1090/ 428-1538/ 792)  24.6 E06-W29 S33-S16
      7 1179  603 (S22W03 L= 78)   762    290        1 (1081/ 531-1264/ 710)  34.7 E04-W09 S25-S16
      8  697  591 (S16E28 L= 46)   548    122        1 ( 669/ 530- 738/ 647)  49.4 E31-E24 S19-S12
      9 1654  539 (S37W41 L=116)   438     84        0 (1619/ 520-1688/ 570)   8.1 W37-W44 S38-S35
     10 1427  452 (S38W20 L= 94)   615    167        0 (1364/ 403-1497/ 507)  24.4 W14-W25 S41-S35
     11  794  443 (S27E26 L= 48)   562     98        1 ( 772/ 398- 821/ 490)  74.1 E28-E23 S30-S23
  

    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