Search found 24 matches

by AntonioG
Mon Aug 20, 2018 3:26 pm
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 03: Filtrado espacial de señales
Replies: 0
Views: 2

Semana 03: Filtrado espacial de señales

Temas a cubrir: Filtrado espacial iltros de media, filtro gausssiano, mediana, , alpha-media, media geométrica, ewma, filtro adaptativo MMSE<-fi-(sn^2/s_L^2)*(f-m) mínimo, máximo, cuantiles, morfologicos varianza, desviacion estándar Filtros de derivadas (primer orden, segundo orden) Filtros lineale...
by AntonioG
Fri Aug 17, 2018 6:53 pm
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 02: Ecualización de histogramas
Replies: 5
Views: 43

Ecualización de imágenes a color

library(jpeg) ruta<-file.choose() img <- readJPEG(ruta) nr<-dim(img)[1] nc<-dim(img)[2] kernelF <- function(x) { m<-length(x) res<-numeric() for(i in 1:m) { r<-integrate(f,0,x[i],stop.on.error = FALSE) res<-c(res,min(r$value,1)) } return(res) } T<-function(x){ return(kernelF(x)) } Ecualiza<-functio...
by AntonioG
Fri Aug 17, 2018 2:34 am
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 02: Ecualización de histogramas
Replies: 5
Views: 43

Kernel density estimation

library(jpeg) ruta<-file.choose() img <- readJPEG(ruta) R <- img[,,1] plot(0:1,0:1, type='n', axes = FALSE, xlab = "", ylab = "") rasterImage(R, 0, 0, 1, 1) kernelf <- function(x) { m<-length(x) h <- 0.00725 res<-numeric() for(i in 1:m) { r <- sum(dnorm((x[i] - xi)/h) / (length(xi)*h)) res<-c(res,r...
by AntonioG
Fri Aug 17, 2018 2:33 am
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 02: Ecualización de histogramas
Replies: 5
Views: 43

Gaussian Mixture Model

# Algorirmo EM para GMM http://www.aishack.in/tutorials/expectation-maximization-gaussian-mixture-model-mixtures/ library(jpeg) ruta<-file.choose() img <- readJPEG(ruta) R <- img[,,1] plot(0:1,0:1, type='n', axes = FALSE, xlab = "", ylab = "") rasterImage(R, 0, 0, 1, 1) library(mixtools) xi<-as.vec...
by AntonioG
Thu Aug 16, 2018 7:36 pm
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 02: Ecualización de histogramas
Replies: 5
Views: 43

F empirica

# EStimacion F empirica library(jpeg) ruta<-file.choose() img <- readJPEG(ruta) R <- img[,,1] plot(0:1,0:1, type='n', axes = FALSE, xlab = "", ylab = "") rasterImage(R, 0, 0, 1, 1) T<-function(x){ # F empirica F<-sum(R<x)/prod(dim(R)) return(F) } nr<-dim(R)[1] nc<-dim(R)[2] Res<-matrix(0, nrow=nr, ...
by AntonioG
Mon Aug 13, 2018 2:53 pm
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 02: Ecualización de histogramas
Replies: 5
Views: 43

Código base

Code: Select all

library(jpeg)

ruta<-file.choose()
img <- readJPEG(ruta)

R <- img[,,1]

plot(0:1,0:1, type='n', axes = FALSE, xlab = "", ylab = "")
rasterImage(R, 0, 0, 1, 1)
by AntonioG
Mon Aug 13, 2018 4:45 am
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 02: Ecualización de histogramas
Replies: 5
Views: 43

Semana 02: Ecualización de histogramas

Temas a cubrir: Mejora de contraste como cambio de paleta de tonos de gris. Algoritmo básico de ecualización de histogramas HE. Método de la transformada inversa Algoritmos y AHE. Otras formas de estimar la distribución de las intensidades (F empirica, estimación de densidad por métodos kernel y mez...
by AntonioG
Fri Aug 10, 2018 6:49 pm
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 01: Formación de señales e imágenes digitales
Replies: 6
Views: 48

Visualización de la imagen digital de ejemplo

ruta<-file.choose() x<-read.table(ruta,header=FALSE,sep=" ") str(x) R<-matrix(x$V1,ncol=8,byrow=TRUE) G<-matrix(x$V2,ncol=8,byrow=TRUE) B<-matrix(x$V3,ncol=8,byrow=TRUE) R<-R[nrow(R):1,] G<-G[nrow(R):1,] B<-B[nrow(R):1,] par(mar=c(0, 0, 0, 0)) image(R, useRaster=TRUE, axes=FALSE) RGB<-array(0,dim=c...
by AntonioG
Fri Aug 10, 2018 6:48 pm
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 01: Formación de señales e imágenes digitales
Replies: 6
Views: 48

Ejemplo de imagen digital de 8x8

1769 2210 764 833 1163 431 919 1315 496 1026 1511 587 974 1380 522 957 1327 486 958 1357 511 1427 2467 908 929 1259 470 1012 1441 551 1131 1706 676 1212 1856 770 1024 1472 578 839 1146 432 894 1255 463 1746 3390 1204 926 1244 463 1597 2606 1165 1468 2350 1016 1618 2627 1155 1115 1507 574 1045 1413 ...
by AntonioG
Wed Aug 08, 2018 3:31 pm
Forum: Procesamiento de señales e imágenes digitales
Topic: Semana 01: Formación de señales e imágenes digitales
Replies: 6
Views: 48

Teorema del muestreo

f<-function(t) { ft<-0.5*sin(25*2*pi*t)+cos(13*2*pi*t) } curve(f,from=0,to=1,n=2000,lwd=2) k<-51 T<-1/k samples<-seq(0,1,T) xn<-f(samples) plot(xn,type="h") points(xn,pch=19) curve(f,from=0,to=1,n=2000,lwd=2) points(samples,xn,pch=19) rect<-function(x) { res<-ifelse(x>=0 & x<1,1,0) return(res) } tr...