Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1373 1631 (N42W34 L=289)   341     59        1 (1341/1616-1403/1647) -31.0 W31-W37 N40-N43
      2 1244 1606 (N41W22 L=277)   295     58        1 (1212/1581-1262/1634)  38.7 W18-W24 N39-N43
      3 1781 1260 (N14W58 L=313)  1423    217        2 (1706/1176-1837/1331) -59.0 W52-W64 N08-N18
      4  315 1164 (N14E51 L=204)   384     85        1 ( 285/1145- 348/1193) -18.4 E53-E47 N13-N16
      5 1096  862 (S06W04 L=259)   335    121        0 (1031/ 830-1143/ 887) -26.6 W00-W06 S08-S04
      6  419  841 (S05E42 L=213)   260     76        1 ( 405/ 816- 439/ 881) -63.4 E42-E40 S06-S02
      7 1436  640 (S22W26 L=281)   562     73        2 (1410/ 603-1472/ 672) -68.2 W24-W29 S24-S19
      8  575  594 (S22E33 L=222)   408     86        1 ( 544/ 563- 611/ 624) -45.0 E35-E31 S23-S19
      9 1306  529 (S29W18 L=273)   726    155        2 (1248/ 482-1364/ 574)  36.9 W13-W21 S32-S26
     10  948  499 (S29E08 L=247)   753    283        0 ( 844/ 439-1066/ 556)  26.6 E16-W00 S33-S25
     11  553  516 (S27E37 L=218)   372    106        0 ( 529/ 471- 593/ 552)  45.0 E40-E33 S30-S24
     12 1576  392 (S43W48 L=303)   647     43        1 (1551/ 372-1606/ 416) -24.4 W43-W54 S45-S40
  

    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