Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  821 1725 (N33E31 L= 94)   863    160        2 ( 790/1651- 851/1789) -39.3 E36-E25 N28-N38
      2  511 1744 (N30E59 L= 66)  1177    140        1 ( 469/1691- 555/1788) -16.3 E73-E49 N26-N33
      3 1236 1615 (N34W01 L=126)   292     72        1 (1211/1593-1269/1642) -18.4 W00-W03 N32-N36
      4  688 1602 (N22E34 L= 91)   295    108        1 ( 645/1569- 739/1639)  51.3 E36-E31 N19-N25
      5  406 1572 (N16E56 L= 69)   598    114        1 ( 378/1542- 443/1615) -18.4 E59-E50 N14-N19
      6  980 1485 (N20E12 L=113)  1348    446        7 ( 853/1346-1115/1608) -28.6 E22-W00 N13-N25
      7  726 1385 (N09E25 L=100)  1363    283        3 ( 700/1219- 758/1579) -70.0 E30-E20 S01-N21
      8 1498 1419 (N27W22 L=147)   498    125        0 (1478/1378-1522/1480)  82.9 W20-W24 N23-N31
      9 1627 1402 (N29W33 L=158)   675    206        0 (1557/1315-1673/1472) -38.7 W26-W37 N24-N32
     10 1121 1397 (N17E02 L=123)   307     56        1 (1089/1378-1155/1417)  33.7 E03-W00 N15-N17
     11 1043  910 (S13W03 L=128)   637    197        1 ( 983/ 884-1143/ 949)  32.0 W00-W08 S15-S10
     12  617  827 (S26E22 L=103)   301     85        1 ( 568/ 800- 651/ 851)  -9.5 E25-E19 S26-S24
     13  568  690 (S35E25 L=100)   336     79        1 ( 539/ 667- 607/ 718)  -9.5 E27-E21 S36-S33
     14 1418  501 (S28W40 L=165)   913    190        3 (1338/ 441-1484/ 572) -10.6 W31-W47 S30-S25
     15 1245  523 (S31W27 L=152)   375     95        0 (1197/ 508-1296/ 537)  15.9 W23-W30 S32-S29
     16  733  486 (S46E09 L=116)   808    144        3 ( 674/ 460- 819/ 518)  18.4 E14-E02 S48-S43
     17 1419  370 (S36W48 L=173)   604     82        1 (1384/ 341-1464/ 399)  68.2 W45-W50 S38-S33
     18  744  234 (S68E03 L=122)   916    147        0 ( 683/ 205- 813/ 264)  -2.6 E13-W09 S70-S65
  

    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