Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  131 1211 (S13E74 L=291)   945    141        0 ( 112/1169- 158/1246)  59.0 E76-E70 S15-S10
      2 1681 1139 (N20W39 L= 43)   691    149        1 (1632/1111-1742/1166)  11.3 W34-W44 N18-N21
      3  322 1055 (S18E47 L=318)   404    126        0 ( 273/1023- 368/1075)  32.0 E50-E42 S21-S16
      4 1524  980 (N07W30 L= 35)   343     77        0 (1492/ 932-1547/1023)  80.5 W29-W31 N03-N09
      5  882  760 (S22E01 L=  3)   518     57        1 ( 846/ 721- 918/ 793)  68.2 E03-W00 S25-S20
      6 1241  709 (S16W21 L= 26)   330    107        0 (1176/ 688-1292/ 734)  39.8 W17-W23 S18-S13
      7 1218  561 (S25W25 L= 30)   322     97        0 (1159/ 547-1264/ 576)  15.9 W21-W28 S26-S23
      8  838  500 (S40W03 L=  8)   916    172        2 ( 788/ 427- 904/ 571) -26.6 E02-W10 S43-S36
      9 1033  432 (S38W20 L= 24)   367     81        0 ( 986/ 412-1076/ 453)   8.1 W16-W23 S39-S37
     10 1533  417 (S24W54 L= 59)   733    210        1 (1460/ 398-1609/ 429)  22.6 W48-W60 S26-S21
  

    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