Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1064 1826 (N48E17 L=195)  1220    275        2 (1027/1744-1105/1896) -74.5 E20-E14 N39-N57
      2  354 1508 (N13E56 L=157)   452    102        0 ( 324/1463- 398/1534)  51.3 E58-E52 N09-N14
      3 1438 1414 (N25W19 L=231)   662    109        0 (1403/1354-1475/1464)  63.4 W17-W21 N20-N28
      4  352 1380 (N05E51 L=161)   579    120        0 ( 329/1317- 388/1434)  82.9 E53-E49 N00-N08
      5 1734 1293 (N24W43 L=255)  1062    132        1 (1690/1221-1767/1364) -63.4 W38-W45 N20-N28
      6  754 1305 (N05E21 L=192)   200     48        0 ( 728/1281- 777/1332)  33.7 E22-E19 N03-N07
      7  782  967 (S14E13 L=199)  1080    309        3 ( 672/ 829- 867/1122)  70.7 E18-E10 S23-S03
      8 1616 1049 (N06W36 L=248)   390    104        0 (1570/1014-1654/1097) -20.6 W31-W39 N04-N08
      9  325  941 (S22E46 L=167)  1105    334        0 ( 225/ 840- 405/1035)  40.9 E54-E39 S28-S15
     10  423  778 (S31E37 L=176)  1414    396        0 ( 299/ 691- 551/ 893)  36.9 E48-E27 S37-S21
     11 1595  653 (S16W45 L=257)  4549    979        4 (1156/ 551-1871/ 744)  20.8 W18-W76 S31-S05
  

    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