Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1270 1760 (N55W29 L=167)   785    123        1 (1214/1728-1332/1789)  10.0 W22-W39 N52-N57
      2  394 1464 (N32E52 L= 86)  1738    331        2 ( 308/1352- 505/1576) -23.6 E70-E37 N24-N38
      3 1147 1549 (N38W10 L=148)   597    142        1 (1089/1517-1216/1573) -18.4 W06-W15 N35-N39
      4 1328 1503 (N34W24 L=161)   329     77        0 (1293/1487-1366/1518)   8.1 W20-W27 N33-N35
      5  846 1246 (N18E11 L=126)   616    193        0 ( 802/1168- 913/1302)  58.0 E14-E06 N12-N21
      6  146 1224 (N14E76 L= 61)  1563     90        1 ( 123/1178- 174/1257)  51.3 E83-E71 N11-N16
      7  549 1024 (N04E31 L=107)  1105    265        1 ( 472/ 933- 653/1100)  36.9 E36-E23 S01-N08
      8  480  870 (S06E36 L=101)   362     45        1 ( 457/ 848- 512/ 892) -14.0 E37-E33 S07-S04
      9  428  624 (S22E44 L= 93)  2311    399        3 ( 332/ 447- 497/ 736)  35.0 E69-E36 S35-S14
     10 1741  700 (S19W54 L=191)   389     75        1 (1703/ 677-1764/ 723)  18.4 W49-W56 S20-S17
     11 1186  491 (S31W11 L=148)  3936    698        9 ( 930/ 342-1419/ 714)  38.2 E08-W25 S43-S16
     12  821  507 (S30E15 L=122)  1626    484        0 ( 673/ 381-1013/ 657)  34.6 E30-E01 S39-S19
  

    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