Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  910 1895 (N52E37 L=292)   900    102        2 ( 872/1869- 955/1921)  -9.5 E44-E31 N48-N54
      2  666 1795 (N37E48 L=281)   878     92        3 ( 636/1750- 694/1849) -36.9 E54-E42 N33-N42
      3  812 1729 (N34E31 L=298)   282     84        1 ( 796/1690- 831/1770) -45.0 E33-E27 N31-N37
      4 1249 1442 (N22W06 L=335)   282     92        1 (1216/1410-1284/1466) -21.8 W03-W08 N20-N23
      5 1054 1374 (N14E04 L=325)   204     68        0 (1021/1351-1091/1395)  53.1 E05-E02 N12-N16
      6 1867 1078 (N14W59 L= 28)   596     72        1 (1837/1053-1891/1112)  45.0 W55-W61 N12-N16
      7  451  948 (S20E36 L=293)   516    114        0 ( 417/ 898- 499/ 987)  50.2 E38-E32 S22-S16
      8 1705  646 (S14W55 L= 24)   403     41        0 (1684/ 619-1724/ 674)  90.0 W53-W57 S15-S12
      9  346  580 (S45E50 L=279)   742     78        1 ( 321/ 546- 370/ 625)  78.7 E54-E46 S47-S41
     10 1684  522 (S21W60 L= 29)  1973    257        7 (1611/ 438-1773/ 597)  69.0 W54-W65 S27-S14
     11 1118  528 (S34W17 L=346)   625    220        0 ( 999/ 509-1213/ 551)  24.8 W09-W22 S37-S31
     12 1575  522 (S24W50 L= 19)   309     54        2 (1556/ 494-1594/ 547)  90.0 W48-W51 S26-S22
     13  993  257 (S57W21 L=350)   991    209        1 ( 882/ 216-1079/ 307)   0.0 W07-W31 S58-S54
  

    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