Genome Browser Photo Gateway

#!/usr/bin/env perl # DO NOT EDIT the /cluster/bin/scripts copy of this file -- # edit ~/kent/src/hg/utils/mkPhotoGateway.pl instead. # $Id: mkPhotoGateway.pl,v 1.2 2007/06/01 19:39:14 galt Exp $ use warnings; use strict; sub usage() { print "usage: ./mkPhotoGateway.pl \n"; print "\tthe resulting page is written to the file.html specified\n"; print "\tCurrent contents of dbDb.hgcentraltest are used to construct\n"; print "\tthe resulting gateway page.\n"; } my $argc=scalar(@ARGV); if ($argc != 1) { usage; die "ERROR: Please supply a single file name for results.\n"; } my $resultFile = shift; open (FH,">$resultFile") or die "ERROR: can not write to $resultFile"; my @dbs = split('\n', `hgsql -N -e "select name from dbDb order by orderKey;" hgcentraltest`); my @sciNames = split('\n', `hgsql -N -e "select scientificName from dbDb order by orderKey;" hgcentraltest`); die "incorrect db vs. sciName count" if (scalar(@dbs) != scalar(@sciNames)); printf STDERR "in dbDb: dbs: %d == %d scientific names\n", scalar(@dbs), scalar(@sciNames); my @dbImages; for (my $i = 0; $i < scalar(@dbs); ++$i) { chomp ($dbs[$i]); chomp ($sciNames[$i]); next if ($dbs[$i] =~ m/hg100/); next if ($dbs[$i] =~ m/canHg/); $sciNames[$i] =~ s/ /_/g; my $jpgName = "$sciNames[$i].jpg"; if ( -s "/usr/local/apache/htdocs/images/$jpgName" ) { push @dbImages,$dbs[$i]; } } my @useDbs; my $prevDb = ""; my $prevVer = -1; my $latestDb = "xxx"; for (my $i = 0; $i < scalar(@dbImages); ++$i) { next if ($dbImages[$i] =~ m/^rheMac[0-9][a-z]$/); next if ($dbImages[$i] =~ m/^pt0$/); if (length($prevDb) > 0) { $dbImages[$i] =~ m/[0-9]+$/; my $dbVersion = $&; my $stripDigits = $`; if ($stripDigits =~ m/^$prevDb$/) { if ($dbVersion > $prevVer) { $prevVer = $dbVersion; $latestDb = $dbImages[$i]; } } else { push @useDbs,$latestDb; $prevVer = $dbVersion; $latestDb = $dbImages[$i]; } $prevDb = $stripDigits; } else { $dbImages[$i] =~ m/[0-9]+$/; $prevVer = $&; $latestDb = $dbImages[$i]; $prevDb = $dbImages[$i]; $prevDb =~ s/[0-9]+$//; } } push @useDbs,$latestDb; printf STDERR "found %d images, %d latest dbs\n", scalar(@dbImages), scalar(@useDbs); print FH "Genome Browser Photo Gateway\n"; print FH '', "\n"; print FH "\n\n"; print FH "\n"; print FH '', "\n"; print FH '', "\n"; my $done = 0; while ($done < scalar(@useDbs)) { print FH "\n"; for (my $i = 0; $i < 3; ++$i) { if ($done < scalar(@useDbs)) { my $href=`hgsql -N -e 'select genome,name from dbDb where name="$useDbs[$done]";' hgcentraltest`; chomp $href; my $genomeOrg = $href; $genomeOrg =~ s/\t/ - /; $href =~ s/\t/&db=/; $href =~ s/ /+/g; my $sciName=`hgsql -N -e 'select scientificName from dbDb where name="$useDbs[$done]";' hgcentraltest`; chomp $sciName; $sciName =~ s/ /_/g; printf " ", $genomeOrg; printf FH " \n", $href; ++$done; } else { printf " "; print FH " \n"; } } printf "\n"; print FH "\n"; } print FH "\n"; close (FH);

- - - - - - - -

\n", $href; printf FH "

", $sciName; printf FH " %s
\n", $genomeOrg; printf FH " (photo credits)