Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1051 1832 (N48E19 L=195)  1110    238        2 (1014/1763-1091/1897) -71.6 E21-E15 N41-N56
      2  341 1385 (N05E53 L=161)   650    157        0 ( 316/1318- 379/1474)  90.0 E53-E51 N00-N10
      3 1421 1419 (N24W18 L=232)   688    138        2 (1387/1356-1460/1473)  71.6 W15-W20 N19-N28
      4 1719 1296 (N24W41 L=255)   905    111        1 (1679/1231-1752/1356) -54.5 W37-W44 N20-N27
      5  703 1229 (S00E22 L=191)   208     91        0 ( 654/1204- 751/1250)  26.6 E25-E19 S02-N01
      6 1598 1058 (N06W34 L=248)   405    112        0 (1552/1023-1639/1110) -23.2 W30-W37 N04-N08
      7  759  962 (S15E15 L=199)   830    284        2 ( 650/ 842- 823/1096)  67.6 E20-E11 S22-S05
      8  317  947 (S22E46 L=168)   837    253        3 ( 246/ 861- 382/1042)  50.2 E52-E41 S27-S15
      9  402  773 (S31E39 L=175)  1186    399        0 ( 290/ 693- 531/ 893)  36.9 E49-E28 S37-S22
     10 1610  667 (S15W45 L=259)  3871    825        7 (1227/ 599-1869/ 726)  18.1 W21-W76 S27-S06
  

    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