Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  694 1713 (N41E29 L=135)   461     74        1 ( 652/1685- 729/1738)  38.7 E31-E26 N38-N43
      2  294 1589 (N32E61 L=103)   843    170        0 ( 268/1561- 327/1608) -18.4 E65-E56 N30-N34
      3 1691 1304 (N16W44 L=208)   280     56        0 (1668/1277-1716/1330) -56.3 W42-W46 N14-N17
      4  180 1014 (S04E61 L=103)   378     93        1 ( 152/ 997- 208/1027)   9.5 E64-E57 S05-S03
      5 1157  871 (S13W08 L=172)   307    112        0 (1098/ 835-1215/ 913)  33.7 W04-W11 S14-S09
      6 1410  844 (S13W25 L=188)   258     80        0 (1389/ 795-1436/ 886)  71.6 W23-W26 S16-S10
      7  181  697 (S23E65 L= 99)  1186    246        3 ( 154/ 650- 220/ 732)  21.8 E69-E59 S25-S21
      8 1857  665 (S21W66 L=230)   727    182        5 (1827/ 615-1896/ 709)  50.2 W62-W70 S24-S18
      9  455  576 (S32E41 L=123)   745    258        2 ( 393/ 530- 520/ 658)  31.6 E47-E34 S35-S26
     10  891  409 (S44E09 L=155)   351     62        2 ( 857/ 391- 929/ 431)  -9.5 E11-E05 S45-S42
     11 1228  363 (S46W20 L=183)   351     62        0 (1202/ 339-1261/ 393)  45.0 W17-W21 S48-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