Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  921 1888 (N64W33 L=358)  1855    230        3 ( 846/1814-1005/1962)   3.9 W19-W65 N57-N67
      2 1196 1740 (N41W37 L=  1)  2293    426        4 (1011/1699-1396/1783) -26.6 W21-W50 N33-N46
      3 1664 1620 (N20W65 L= 29)   813     36        1 (1631/1596-1693/1649) -31.0 W60-W69 N18-N21
      4  799 1592 (N43W01 L=325)   374     95        0 ( 777/1542- 825/1631)  56.3 E01-W04 N39-N45
      5 1751 1493 (N11W64 L= 28)   629     89        1 (1721/1444-1775/1533)  80.5 W62-W66 N07-N13
      6 1572 1109 (S05W34 L=359)  2152    453        2 (1519/ 888-1606/1310)  60.3 W26-W40 S15-N06
      7  314 1097 (N23E45 L=279)  1877    294        6 ( 214/1008- 402/1187) -39.1 E55-E38 N16-N30
      8 1884 1154 (S11W63 L= 28)   432     51        0 (1858/1134-1911/1175) -26.6 W60-W66 S12-S09
      9  133 1110 (N28E68 L=257)  1784    131        2 (  83/1051- 180/1160) -22.8 E81-E60 N23-N31
     10  744 1085 (N14E15 L=310)   794    228        0 ( 656/1032- 839/1160) -54.0 E18-E09 N08-N19
     11  270  948 (N16E51 L=273)   324     98        1 ( 242/ 907- 313/ 979)  26.6 E54-E46 N13-N16
     12  336  867 (N10E47 L=278)   337     62        0 ( 307/ 833- 356/ 888) -45.0 E49-E45 N07-N10
     13  226  800 (N08E60 L=265)   469     91        1 ( 207/ 767- 257/ 845) -68.2 E61-E56 N05-N10
     14  298  793 (N06E52 L=272)   334     44        0 ( 278/ 775- 324/ 816) -63.4 E54-E49 N05-N07
     15  252  725 (N03E59 L=265)   893    184        3 ( 193/ 675- 289/ 776) -54.5 E64-E55 N00-N07
     16 1265  606 (S26W04 L=329)   570    195        1 (1206/ 539-1333/ 669)  26.6 E01-W09 S28-S23
     17 1751  570 (S44W45 L= 10)  1094    154        2 (1686/ 524-1819/ 623)   7.1 W37-W54 S45-S40
     18  737  256 (S36E46 L=279)   852    128        0 ( 688/ 216- 782/ 316)  17.4 E53-E37 S38-S32
  

    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