Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  920 1695 (N30E28 L=326)   361     97        1 ( 865/1679- 972/1709)  33.7 E30-E24 N27-N31
      2 1470 1628 (N45W11 L=  5)   797    116        1 (1428/1592-1507/1657) -26.6 W07-W14 N41-N47
      3 1540 1606 (N46W18 L= 12)   470    108        0 (1525/1561-1562/1639)  82.9 W16-W19 N41-N49
      4 1690 1428 (N38W33 L= 27)  2960    378        9 (1581/1279-1765/1543) -47.1 W22-W41 N29-N44
      5  390 1484 (N04E55 L=299)   513    126        0 ( 354/1446- 414/1510) -45.0 E58-E51 N01-N05
      6 1041 1259 (N07E05 L=348)   283     90        1 (1026/1219-1063/1309) -59.0 E06-E03 N04-N09
      7 1754 1248 (N29W41 L= 35)   723    113        3 (1724/1209-1777/1297)  74.1 W38-W42 N25-N32
      8  172 1266 (S12E68 L=286)   941    193        1 ( 132/1225- 209/1295)  11.3 E74-E64 S14-S10
      9 1553 1154 (N16W27 L= 21)   337     92        0 (1519/1105-1571/1202) -33.7 W23-W29 N13-N17
     10 1803 1102 (N22W48 L= 41)   509    114        0 (1770/1079-1839/1135)   8.1 W44-W51 N20-N23
     11 1464 1016 (N05W25 L= 19)   882    217        2 (1395/ 964-1554/1099)  -8.7 W18-W31 N03-N07
     12 1868  981 (N17W57 L= 51)   962    124        2 (1847/ 923-1886/1042) -45.0 W54-W61 N14-N20
     13 1710  942 (N09W44 L= 38)  1245    258        2 (1615/ 912-1794/ 968)  26.6 W37-W51 N04-N12
     14 1720  760 (N00W51 L= 45)  1049    174        0 (1649/ 699-1774/ 840)  63.4 W47-W53 S05-N05
     15 1258  792 (S13W20 L= 14)   520    200        1 (1158/ 771-1337/ 839)  37.9 W14-W23 S16-S08
     16  731  753 (S30E10 L=344)   212     55        0 ( 700/ 734- 755/ 780)  63.4 E11-E08 S31-S27
     17 1132  639 (S25W18 L= 12)  1126    296        3 ( 996/ 536-1257/ 715)  -5.2 W07-W29 S27-S21
     18 1523  597 (S15W44 L= 37)   492    129        0 (1471/ 557-1584/ 632)  58.0 W41-W46 S18-S10
  

    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