Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  680 1748 (N29E50 L=235)  1494    143        2 ( 637/1703- 735/1817) -21.3 E59-E41 N25-N34
      2  395 1595 (N13E63 L=222)  3241    363        7 ( 325/1417- 462/1716) -58.2 E74-E56 N00-N21
      3 1701 1562 (N53W41 L=326)  1081    138        1 (1659/1506-1749/1596) -18.4 W35-W47 N49-N56
      4  555 1539 (N13E45 L=240)   377     87        0 ( 533/1499- 585/1577) -78.7 E46-E43 N10-N15
      5 1583 1281 (N26W28 L=312)   244     82        0 (1554/1245-1607/1318) -26.6 W24-W30 N24-N27
      6 1738 1216 (N27W42 L=327)   465     69        1 (1714/1171-1761/1246) -53.1 W39-W44 N24-N29
      7  285 1091 (S19E52 L=233)   340     55        1 ( 263/1062- 306/1113) -21.8 E53-E48 S20-S17
      8 1751  916 (N10W49 L=334)   415    100        0 (1739/ 870-1764/ 969) -59.0 W47-W52 N08-N13
      9  418  812 (S33E35 L=249)   541    146        0 ( 361/ 777- 480/ 849)  38.7 E40-E30 S36-S28
     10  826  602 (S34W00 L=285)   361     59        2 ( 793/ 579- 863/ 623)  53.1 E01-W02 S36-S32
     11 1586  496 (S17W55 L=340)  3384    412        2 (1497/ 333-1734/ 607) -79.2 W50-W66 S27-S06
     12 1330  423 (S29W41 L=326)   322     78        0 (1284/ 411-1368/ 435)  38.7 W38-W43 S31-S27
  

    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