Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  548 1610 (N23E49 L= 41)   472    116        1 ( 502/1586- 609/1633)   0.0 E53-E42 N22-N24
      2 1617 1457 (N39W34 L=125)   428    103        0 (1577/1429-1658/1481)  -7.1 W30-W38 N37-N40
      3 1334 1398 (N27W11 L=102)   233    152        0 (1281/1370-1394/1413)  14.0 W08-W16 N24-N28
      4 1222 1362 (N22W05 L= 95)   244     43        2 (1200/1338-1244/1389) -33.7 W03-W06 N20-N24
      5 1086 1331 (N17E03 L= 87)   727    244        4 ( 995/1263-1156/1381)   6.3 E08-W01 N14-N19
      6 1325 1304 (N21W12 L=103)   353     47        0 (1297/1280-1349/1332)  45.0 W10-W14 N19-N23
      7 1332 1231 (N17W14 L=104)   539    136        1 (1283/1197-1375/1256)  -8.1 W10-W17 N14-N18
      8  961  587 (S29W06 L= 96)  1309    376        1 ( 790/ 530-1113/ 653)  37.9 E04-W14 S36-S21
      9  834  357 (S50W04 L= 95)   425     58        0 ( 805/ 337- 879/ 379)  15.9 W01-W08 S51-S48
     10 1108  289 (S47W31 L=121)  1752    229        2 (1000/ 263-1216/ 336)  22.4 W21-W40 S50-S43
  

    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