Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1633 1506 (N38W45 L=131)  1383    256        2 (1540/1402-1702/1593) -50.2 W37-W49 N30-N44
      2  948 1507 (N32E10 L= 76)   313     55        1 ( 922/1489- 978/1535) -36.9 E11-E07 N30-N33
      3 1095 1403 (N26W01 L= 87)   216     41        1 (1071/1383-1118/1419)   0.0 W00-W03 N24-N27
      4  203 1319 (N12E68 L= 18)  1004    149        3 ( 178/1288- 233/1355) -53.1 E72-E62 N10-N14
      5  838  808 (S13E10 L= 76)   333    102        1 ( 819/ 747- 857/ 856) -63.4 E11-E08 S16-S10
      6  693  778 (S16E19 L= 66)   255    149        0 ( 659/ 728- 753/ 805)  33.7 E21-E15 S19-S14
      7 1641  718 (S13W45 L=131)   327     91        1 (1616/ 687-1679/ 738)  45.0 W43-W48 S14-S11
      8 1184  676 (S19W13 L= 99)   219     86        1 (1160/ 645-1212/ 715)  78.7 W11-W14 S20-S15
      9 1715  675 (S15W53 L=139)   482    124        3 (1678/ 638-1752/ 701)  38.7 W49-W56 S17-S13
     10  431  647 (S27E40 L= 46)   467    127        0 ( 389/ 618- 483/ 696) -33.7 E43-E35 S29-S24
     11 1716  579 (S21W57 L=142)   811    204        2 (1644/ 549-1764/ 622)  20.0 W50-W62 S23-S18
     12 1311  315 (S44W33 L=119)   361     83        1 (1278/ 295-1342/ 336)  -8.1 W29-W36 S45-S42
  

    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