Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1362 1573 (N41W13 L=130)   261     57        1 (1332/1556-1395/1588)   9.5 W10-W16 N39-N41
      2 1563 1493 (N39W31 L=148)   358     96        0 (1524/1476-1601/1507)  15.9 W27-W34 N37-N40
      3 1259 1384 (N25W08 L=125)   698    138        1 (1204/1329-1331/1445) -21.8 W03-W13 N22-N28
      4 1776 1394 (N37W51 L=168)  1157    195        1 (1724/1349-1828/1427) -14.0 W45-W56 N34-N39
      5 1720 1071 (N15W46 L=163)  1403    306        3 (1671/ 948-1787/1153)  85.6 W42-W51 N07-N20
      6 1463  824 (S05W31 L=148)   439    116        0 (1410/ 806-1513/ 852)   0.0 W26-W34 S05-S03
      7  205  775 (S31E58 L= 59)   629     93        0 ( 179/ 749- 241/ 805)  -8.1 E60-E53 S32-S29
      8 1493  561 (S20W41 L=158)   339     88        1 (1462/ 535-1520/ 587)  -9.5 W37-W43 S21-S18
  

    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