Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1470 1672 (N48W12 L= 77)   316     31        1 (1445/1654-1494/1689)   0.0 W09-W14 N46-N49
      2 1545 1550 (N40W20 L= 86)   752    123        2 (1506/1499-1578/1609) -76.0 W16-W23 N36-N44
      3  244 1463 (N02E70 L=356)   757     56        1 ( 214/1439- 282/1482)  -4.4 E77-E64 N00-N02
      4 1027 1418 (N15E10 L= 55)   318    103        0 ( 998/1365-1046/1463) -39.8 E13-E07 N12-N17
      5 1239 1218 (N09W07 L= 72)   187     78        0 (1227/1176-1252/1257) -53.1 W05-W08 N07-N11
      6 1648 1210 (N21W33 L= 99)   357     57        1 (1620/1181-1671/1240) -56.3 W31-W35 N19-N22
      7 1906 1139 (N27W62 L=128)  1750    228        0 (1875/1022-1929/1238) -77.0 W57-W67 N19-N32
      8  219 1135 (S18E59 L=  6)  1159    208        0 ( 162/1099- 301/1182)  14.9 E67-E51 S20-S15
      9 1844 1107 (N22W53 L=119)   474     84        1 (1829/1062-1864/1151) -59.0 W51-W56 N19-N24
     10 1355 1046 (N03W18 L= 83)   262    106        0 (1313/1017-1390/1082)  59.0 W15-W19 N00-N05
     11  187  876 (S35E62 L=  3)   776     98        0 ( 150/ 845- 235/ 913)  -8.7 E68-E55 S36-S32
  

    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