Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  791 1712 (N33E38 L=269)  1238    207        2 ( 726/1645- 866/1786) -20.2 E48-E29 N30-N37
      2  568 1634 (N22E49 L=257)   895    189        0 ( 521/1557- 594/1692) -33.7 E57-E44 N17-N26
      3 1397 1632 (N46W10 L=316)   745    143        0 (1350/1593-1455/1660)  -5.7 W05-W15 N43-N48
      4  440 1570 (N16E56 L=250)   449     84        1 ( 420/1542- 467/1620) -45.0 E61-E52 N13-N18
      5  623 1546 (N18E40 L=267)   465    103        0 ( 588/1501- 651/1593)  90.0 E41-E38 N14-N21
      6 1773 1196 (N27W48 L=355)  1010    127        2 (1743/1128-1801/1254) -49.4 W44-W51 N22-N29
      7 1317 1171 (N12W13 L=320)   274     90        0 (1292/1126-1339/1217) -45.0 W11-W15 N10-N14
      8 1860 1149 (N26W60 L=  6)   703    112        0 (1823/1106-1895/1195) -20.0 W54-W65 N24-N28
      9  563 1077 (S11E29 L=278)   238     54        0 ( 538/1060- 594/1097)   0.0 E31-E26 S11-S09
     10  443  748 (S34E33 L=273)   282     60        1 ( 423/ 732- 478/ 773) -11.3 E35-E30 S34-S32
     11 1281  581 (S22W28 L=335)   325    123        0 (1266/ 532-1308/ 640) -45.0 W25-W31 S24-S18
     12 1550  546 (S17W48 L=355)   334     62        2 (1528/ 517-1570/ 569) -56.3 W46-W50 S18-S15
  

    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