Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1230 1740 (N54W10 L= 86)   580    103        0 (1177/1709-1278/1764) -10.3 W03-W14 N51-N55
      2  488 1507 (N26E46 L= 31)  1273    234        2 ( 439/1430- 524/1620) -69.0 E52-E43 N20-N33
      3 1329 1443 (N31W17 L= 94)  2880    598        4 (1126/1235-1488/1562) -33.7 W01-W28 N18-N38
      4  510 1405 (N20E40 L= 37)   324     46        2 ( 486/1379- 534/1426) -21.8 E42-E37 N18-N20
      5 1178 1224 (N15W08 L= 85)   361    158        0 (1146/1174-1228/1289) -50.2 W05-W11 N12-N18
      6  472  720 (S23E36 L= 41)   282     72        1 ( 441/ 692- 501/ 752)  38.7 E38-E33 S25-S20
      7 1295  671 (S18W21 L= 98)   625    165        0 (1225/ 621-1370/ 726)  45.0 W17-W25 S21-S13
      8 1712  654 (S15W55 L=132)   575     57        0 (1688/ 633-1739/ 678) -18.4 W52-W58 S16-S13
      9 1647  577 (S21W52 L=129)   449     80        0 (1625/ 549-1676/ 620)  78.7 W50-W53 S22-S17
     10 1708  592 (S19W58 L=135)   426     86        0 (1679/ 565-1754/ 614)  23.2 W54-W62 S21-S17
     11  483  516 (S38E38 L= 39)   315     44        1 ( 461/ 495- 506/ 539)  26.6 E41-E35 S39-S35
     12  419  486 (S41E47 L= 30)   694     57        2 ( 394/ 461- 455/ 519) -26.6 E51-E43 S43-S38
  

    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