Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  678 1628 (N38E35 L=339)  1505    285        3 ( 545/1587- 798/1666)  -4.2 E50-E23 N36-N40
      2  306 1525 (N27E69 L=305)   577     56        0 ( 280/1511- 341/1539)  12.5 E75-E63 N25-N28
      3 1292 1396 (N27W15 L= 29)  1338    359        1 (1212/1298-1348/1508) -53.1 W09-W19 N20-N34
      4  918 1465 (N29E11 L=  3)   238     37        1 ( 895/1448- 942/1489) -33.7 E12-E09 N27-N30
      5  463 1362 (N18E43 L=331)   537    170        0 ( 408/1289- 508/1406) -26.6 E48-E36 N14-N20
      6 1741 1254 (N21W52 L= 66)   375    111        0 (1729/1203-1753/1308)  74.1 W50-W53 N17-N24
      7 1585 1273 (N21W37 L= 51)   312     71        1 (1564/1243-1613/1305) -63.4 W35-W39 N19-N23
      8  929 1122 (N07E07 L=  7)   771    266        3 ( 847/1050-1039/1185)  39.3 E11-W00 N02-N11
      9 1290  564 (S25W22 L= 36)   391    126        0 (1228/ 546-1341/ 585)  -6.3 W17-W26 S26-S24
  

    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