Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1191 1749 (N54W05 L= 84)   345     65        1 (1159/1730-1227/1765)   0.0 W01-W09 N52-N55
      2  474 1519 (N26E48 L= 30)  1378    232        0 ( 418/1438- 508/1632) -60.3 E55-E45 N20-N34
      3 1334 1411 (N29W17 L= 96)  1905    523        2 (1136/1248-1448/1544) -36.5 W02-W25 N19-N36
      4 1765 1415 (N33W59 L=138)   772    115        1 (1717/1364-1802/1461) -32.0 W54-W64 N30-N36
      5 1155 1225 (N15W06 L= 85)   406    166        0 (1120/1175-1206/1294) -50.2 W03-W09 N12-N19
      6 1321 1271 (N19W17 L= 95)   210     81        1 (1293/1252-1355/1289) -11.3 W14-W19 N18-N20
      7 1275  677 (S18W20 L= 98)   586    166        0 (1205/ 628-1346/ 727)  48.8 W15-W23 S21-S13
      8  386  698 (S26E44 L= 35)   290     54        0 ( 364/ 676- 416/ 722) -26.6 E45-E40 S27-S24
      9 1673  621 (S18W52 L=131)  1029    172        1 (1611/ 550-1727/ 677)  60.9 W48-W57 S22-S13
     10 1690  573 (S21W57 L=135)   858    147        1 (1662/ 515-1743/ 613)  74.1 W52-W61 S24-S17
     11  426  500 (S40E45 L= 33)  1023    136        2 ( 377/ 466- 488/ 538)  18.4 E52-E37 S43-S36
  

    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