Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1148 1757 (N54W01 L= 83)   566    118        2 (1089/1733-1205/1776)   0.0 E05-W06 N52-N55
      2  420 1494 (N24E52 L= 30)  2211    383        3 ( 348/1398- 474/1632) -49.1 E59-E44 N18-N33
      3 1286 1412 (N28W14 L= 96)  1834    543        2 (1087/1240-1417/1553) -36.9 E01-W23 N18-N36
      4 1765 1386 (N31W58 L=140)   517     82        0 (1743/1340-1791/1420) -59.0 W55-W60 N28-N33
      5 1095 1237 (N15W02 L= 85)   318    157        0 (1062/1192-1143/1302) -56.3 W00-W05 N12-N18
      6 1214  686 (S18W16 L= 98)   740    194        1 (1141/ 631-1291/ 741)  45.0 W11-W19 S21-S13
      7 1699  654 (S15W54 L=136)   342     65        1 (1672/ 635-1721/ 679)  56.3 W51-W56 S16-S13
      8 1591  581 (S21W46 L=129)   442     72        0 (1569/ 550-1612/ 618)  90.0 W44-W48 S23-S18
      9 1651  564 (S22W53 L=135)   611     98        1 (1629/ 521-1681/ 602) -45.0 W49-W57 S24-S19
     10  432  535 (S37E43 L= 39)   363     66        0 ( 409/ 513- 469/ 556)  26.6 E46-E38 S39-S35
     11  374  502 (S41E53 L= 30)   622     65        2 ( 347/ 484- 405/ 528) -26.6 E57-E47 S42-S38
  

    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