Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1106 1599 (N40W02 L=299)  6042   1141       10 ( 629/1485-1592/1749)  -6.5 E46-W43 N32-N48
      2 1707 1420 (N31W52 L=349)  1599    192        4 (1630/1362-1784/1470) -23.2 W45-W60 N27-N33
      3  423 1367 (N19E45 L=252)   380     73        1 ( 403/1332- 446/1403) -38.7 E47-E42 N16-N21
      4 1542 1313 (N23W34 L=331)   492    156        1 (1523/1248-1562/1400)  84.3 W32-W35 N18-N28
      5  709 1375 (N22E23 L=274)   274     55        1 ( 688/1355- 741/1398) -26.6 E25-E21 N20-N23
      6  167 1203 (N06E69 L=228)  2697    346        3 ( 138/1055- 195/1363) -74.5 E75-E64 S02-N16
      7 1458 1304 (N22W27 L=324)   249     80        0 (1438/1267-1483/1341) -59.0 W26-W29 N19-N24
      8  694 1150 (N07E21 L=276)  1358    369        1 ( 672/ 988- 732/1336)  90.0 E23-E19 S02-N19
      9  264 1229 (N09E57 L=240)   901    147        1 ( 243/1173- 292/1296) -82.9 E59-E52 N05-N13
     10 1861 1236 (N19W66 L=  3)   462     66        0 (1823/1221-1898/1250)   4.4 W60-W73 N18-N20
     11  309  796 (S17E50 L=247)  1235    328        2 ( 282/ 648- 337/ 946)  70.7 E55-E46 S27-S07
     12 1885  917 (S00W67 L=  4)   978     98        2 (1841/ 891-1930/ 943)   8.1 W61-W75 S01-N01
     13  925  707 (S18E04 L=293)   257    120        0 ( 872/ 694- 981/ 730) -15.9 E07-W00 S18-S16
     14 1484  539 (S26W37 L=334)  4481    778       12 (1120/ 405-1762/ 672)  24.5 W12-W59 S37-S15
     15  688  635 (S25E20 L=277)   432    107        0 ( 641/ 606- 744/ 663)  29.7 E23-E16 S26-S22
     16  471  511 (S35E41 L=256)   457     76        1 ( 424/ 497- 508/ 530)  12.5 E45-E36 S36-S33
  

    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