Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1104 1611 (N35W21 L=177)  1292    240        3 (1006/1557-1223/1645) -18.4 W13-W28 N31-N37
      2 1746 1375 (N06W57 L=213)   976    113        2 (1698/1326-1790/1427) -74.1 W53-W60 N02-N09
      3  233 1232 (N31E52 L=104)   718    208        2 ( 178/1209- 278/1268) -26.6 E57-E47 N29-N34
      4 1248  920 (S08W10 L=166)  1398    429        4 (1095/ 720-1344/1121)  37.3 E02-W20 S16-N01
      5 1762 1077 (S11W49 L=204)   295     80        1 (1727/1053-1788/1092) -18.4 W45-W51 S12-S09
      6 1268  700 (S21W07 L=163)   207     43        1 (1236/ 685-1295/ 719)   0.0 W05-W09 S22-S20
      7 1105  650 (S21E04 L=152)   194     73        0 (1073/ 629-1143/ 684) -59.0 E04-E01 S22-S17
      8  895  435 (S29E23 L=133)   747    180        1 ( 874/ 361- 918/ 515)  61.4 E26-E20 S34-S23
      9 1308  456 (S39W05 L=161)   420     93        0 (1263/ 423-1339/ 498)  23.2 W00-W07 S40-S36
     10  512  384 (S24E52 L=104)  1589    320        4 ( 469/ 297- 556/ 493)  63.4 E58-E46 S30-S16
  

    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